Black Hole Remnants in the Early Universe 
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We consider the production of primordial micro black holes (MBH) remnants in the early universe. 
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we find that our model is able to explain the quadrupole anomaly of the CMB power spectrum. 
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I. INTRODUCTION 

Inflation is without doubt the best model to explain 
the observed spatially flat and homogeneous Universe. 
Nevertheless, despite the great successes of the standard 
ACDM model in explaining almost all the data on CMB 
anisotropy as most recently measured by WMAP obser- 
vations, the suppression of the I = 2 quadrupole mode 
still remains a puzzle in the framework of the standard 
ACDM model (for a review on this subject, see e.g. [l|). 

Recently, several authors @, Q have been able to shed 
some light on this region of the CMB power spectrum, by 
investigating the possibility of a pre-inflationary epoch, 
dominated by radiation, instead of the usual inflationary 
vacuum. They found that a pre inflation radiation era 
can produce a suppression of the low k modes of the 
primordial power spectrum, and this in turn affects the 
low / modes of CMB anisotropy power spectrum. In fact, 
although inflation has the effect of washing out the initial 
conditions of the Universe, it happens that, if the present 
Universe is just comparable to the size of the inflated 
region, a pre- inflation era may leave imprints on the CMB 
power spectrum. 

However, these early attempts suffered of an arbitrary 
initial condition in the pre-inflationary era. Also, the 
space of the numerical parameters encoding the initial 
radiation density was merely explored, without stating 
precise criteria for the choice of specific numerical values. 

In the present paper we propose a pre-inflationary 
scenario that is based on the generic micro black hole 
(MBH) production and a minimal set of first princi- 
ples, namely the generalized uncertainty principle (GUP) 
and the holographic principle (HP), that can give rise to 
the suppression of the CMB quadrupole self-consistently 
without the need of arbitrary inputs. Specifically, we con- 
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sider the possibility of production of micro black holes 
in the early pre-inflationary Universe, due to quantum 
fluctuations of the metric field 0, as the seeds for 
the suppression of the inflaton fluctuations. There are 
two salient features of this MBH nucleation. One is 
that the production rate per unit volume of space and 
time is very high at the Planck temperature. To pre- 
vent unphysical over-production of MBH, we invoke the 
holographic principle (HP) to constrain the initial condi- 
tion of MBH production. The other is that the rate of 
such MBH production is a strong function of the back- 
ground temperature. In particular, the rate is exponen- 
tially suppressed when the temperature of the universe 
is sufficiently below the Planck temperature. Inflation 
is in general assumed to start when the temperature of 
the universe reaches the scale of the GUT energy, about 
10 15 - 10 16 GeV. Therefore one expects that the MBH 
production activity would cease long before the onset of 
the inflation, and the MBH would have been totally evap- 
orated and the universe would turn into radiation era be- 
fore the inflation begins. However, when the Generalized 
Uncertainty Principle (GUP) is taken into consideration, 
the complete decay of the nucleated MBH into radia- 
tion is prevented, and we have massive, but inert black 
hole remnants @ populating the pre-inflationary phase 
of the Universe. Furthermore, the nucleation of MBH is 
so efficient and fast that the Universe is put into a mat- 
ter dominated era within a few Planck times, just about 
10 3 t p after the Big Bang (i.e., well before inflation) and 
there it stays until the onset of inflation. Such a pre- 
inflation matter-dominated universe then suppresses the 
initial inflaton fluctuations at the onset of the inflation. 

Accurate numerical simulations allowed us to single 
out almost unique numerical values for the relevant ra- 
diation and matter parameters. We have computed the 
effects of a pre-inflationary matter epoch on the primor- 
dial power spectrum of the quantum fluctuations of a 
scalar field, both analytically and numerically. Our ana- 
lytical solution, also a new feature of the present attempt 
with respect to the previous all-numerical investigations, 
has served as a guide for the more precise numerical com- 
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putations. We have considered three alternative scenar- 
ios. The main model presented in this paper attempts 
to explain the suppression of the quadrupole moment of 
the CMB with a pre-inflationary matter era. In order 
to isolate the cause of the CMB quadrupole anomaly, we 
further examine two variations of this model, one with- 
out the GUP, where the black holes decay into radiation 
completely, and one without any black hole nucleated at 
all, both models resulting into a radiation-dominated era 
before inflation. 

In all three cases the primordial power spectra have 
been fed to the CMBFAST code in order to obtain the 
CMB power spectra, and then compared with each other 
and tested against the last WMAP observational data. 
Boundary conditions have been set in the fully inflation- 
ary epoch, and in so doing we avoided any arbitrary as- 
sumption on the state of radiation in the pre-inflationary 
era. The pre-inflation matter model seems to be the only 
one, among those studied, which is able to describe the 
I = 2 mode suppression, although the radiation model 
still presents a better fitting of the data at high I values. 
This conclusion is widely discussed in the last section of 
the paper, where we also suggest avenues for future re- 
search. We find it remarkable that, based on our ab initio 
model without arbitrary input parameters, our resulting 
suppression of the CMB spectrum agrees well with ob- 
servations. 

The paper is organized as follows. In section II we 
will provide a brief overview of the black hole physics 
under the Generalized Uncertainty Principle. The GUP 
will lead to a new mass-temperature relation and define a 
minimum mass and maximum temperature for the black 
holes. Section III will set up the basic equations gov- 
erning the scenario, on the one hand the absorption and 
emission processes which determine the black hole mass, 
and on the other hand the evolution equations of the 
universe depending on its constitution. We will end up 
with a system of four equations, containing black hole 
mass, black hole density, radiation density and the scale 
factor as variables. At the end of that section we de- 
rive a condition for a pre inflation matter era, and we 
present inflationary solutions. In section IV we will nu- 
merically calculate the black hole production and fix the 
parameters in the evolution equations. Section V will 
deal with the equations needed to be solved to obtain 
the primordial power spectrum of the quantum fluctua- 
tions of a scalar field. First we state some approximate 
analytical solutions obtained via the WKB method, and 
then we present the numerical result of the equations. In 
section VI finally we present the CMB power spectrum of 
the temperature anisotropics obtained by our model, and 
compare it to the two alternative cases with a radiation- 
dominated era before inflation. The conclusions of our 
work are contained in section VII. 

Throughout the paper the Planck length is defined as 
£ 2 = Gh/c 3 , the Planck energy as £ p £ p = he/ '2, and the 
Planck mass as M p = £ p /c 2 . 



II. BLACK HOLE PHYSICS 
A. Generalized Uncertainty Principle 

As it is well known from the classical argument of the 
Heisenberg microscope Q, the size Sx of the smallest 
detail of an object, theoretically detectable with a beam 
of photons of energy E, is roughly given by 



Sx 



he 
2E 



(1) 



since larger and larger energies are required to explore 
smaller and smaller details. 

The research on viable generalizations of the Heisen- 
berg uncertainty principle traces back to many decades 
(see for early approaches M: see for a review @ , and for 
more recent approaches |10l4l3 |). In the last 20 years, 
there have been important studies in string theory [l3j 
suggesting that, in gedanken experiments on high energy 
scattering with high momentum transfer, the uncertainty 
relation should be written as 



(2) 



where l p is the Planck length, and f3£ 2 ~ , where X s is 
the characteristic string length. Since in our high energy 
scattering E ~ cp, the stringy Generalized Uncertainty 
Principle (GUP) can be also written as 



Sx > 



he 

2E 



E 

he ' 



(3) 



where E is the energy of the colliding beams. 
A similar modification of the uncertainty principle has 
been proposed [ll|, [TJJ , on the ground of gedanken scat- 
tering experiments involving the formation of micro black 
holes with a gravitational radius of Rs ~ E. It reads 



Sx > 



|§ for E<£ p 



PR S {E) for E>£ pi 



(4) 



where Rs is the Schwarzschild radius associated with the 
energy E, namely Rs = £ p E/£ p . 

Combining linearly the above inequalities we get 



5x>^+(3Rs(E). 



(5) 



Thus, the GUP originating from micro black hole 
gedanken experiments (MBH GUP) can be written as 



(6) 



Also the stringy inspired GUP (ST GUP, eq.©), using 
the relation £ p l p = hc/2, can be written as 



Sx > 



he 
2E 



0l p - 



E 



(7) 
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where /3 is the deformation parameter, generally believed 
to be of 0(1)- Thus, in 4 dimensions the two principles 
coincide. In 4 + n dimensions, however, they lead to 
remarkably different predictions (see (l5j). 



B. From the uncertainty principle to the 
mass-temperature relation 

Naturally, a modification of the uncertainty relation, 
i.e. of the basic commutators, has deep consequences 
on the quantum mechanics, and on the quantum field 
theory built upon it. The general implementation 
of such commutation rules, as regards Hilbert space 
representation, ultraviolet regularization, or modified 
dispersion relations, has been discussed in a vast 
amount of literature (see [IB] for an incomplete list). 
In the present section, we want to focus on the use 
of (generalized) uncertainty relations to compute the 
basic feature of the Hawking effect, namely the formula 
linking the temperature of the black hole to its mass M. 
The seminal results of Hawking and Unruh [l7|, EH are 
rigorously computed using QFT, based on Heisenberg 
uncertainty principle, on curved space-time. However, 
it has been shown [(| [l9|, [2(| that the full calculation of 
QFT in curved space-time (with standard commutators 
for the ordinary uncertainty principle, or with deformed 
commutators for the GUP) can be safely replaced by 
a computation employing only the (generalized) uncer- 
tainty relation and some basic physical considerations, 
in order to obtain the mass-temperature formula. 

The GUP version of the standard Heisenberg formula 
© is 



dx 



he 

2E 



E 

y 



(8) 



which links the (average) wavelength of a photon to its 
energy E. Conversely, with the relation © one can com- 
pute the energy E of a photon with a given (average) 
wavelength A ~ Sx. 

Following loosely the arguments of Refs. d, [H, fl9M22| , 
we can consider an ensemble of unpolarized photons of 
Hawking radiation just outside the event horizon. From 
a geometrical point of view, it's easy to see that the po- 
sition uncertainty of such photons is of the order of the 
Schwarzschild radius Rs of the hole. An equivalent ar- 
gument comes from considering the average wavelength 
of the Hawking radiation, which is of the order of the ge- 
ometrical size of the hole. By recalling that Rs = £ p m, 
where m = M/M p is the black hole mass in Planck units 
(M p = Ep/c 2 ), we can estimate the photon positional 
uncertainty as 

Sx ~ 2^R S = 2^l p m. (9) 
The proportionality constant \x is of order unity and will 



be fixed soon. With © we can rephrase Eq. 
2/im 



E £ v 



(10) 



According to the equipartition principle the average en- 
ergy E of unpolarized photons of the Hawking radiation 
is linked with their temperature T as 



E 



knT. 



(11) 



In order to fix we consider the semiclassical limit /3 — > 
0, and require that formula (1101) predicts the standard 
semiclassical Hawking temperature: 



Ti 



he 3 



he 



H 



8nGk B M AirksRs 



(12) 



This fixes fi = 7T. 

Defining the Planck temperature T p so that £ p — ksTp/2 
and measuring all temperatures in Planck units as G = 
T/Tp, we can finally cast formula (flO|) in the form 



2m 



27t6 



C2nQ, 



(13) 



where we have defined the deformation parameter £ = 
/3/tt 2 . 

As already mentioned, in the semiclassical limit both 
j3 and C tend to zero and © reduces to the ordinary 
Heisenberg uncertainty principle. In this case Eq. (|13[) 
boils down to 



1 



AirO ' 



(14) 



which is the dimensionless version of Hawking's formula 
CCD- 

As we have seen, a computation of the mass- 
temperature relation for black holes based on the GUP 
has resulted in a modification of the Hawking formula 
for high temperatures. In the next subsection, we sum- 
marize as this leads also to the remarkable prediction of 
black hole remnants (see [6j). 



C. Minimum masses, maximum temperatures 

The standard Hawking formula predicts a complete 
evaporation of a black hole, from an initial mass M down 
to zero mass. As we have seen this is a direct conse- 
quence of the Heisenberg principle. However, when the 
mass-temperature relation is derived from the GUP in- 
stead, the formulation immediately leads to a minimum 
mass and a maximum temperature for the evaporating 
black hole. Precisely we have, for the GUP, 



i 



2nVC 



(15a) 



(15b) 
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Note that, as expected, max — > °o and m min — > in 
the Hawking limit j3 — > 0. Therefore the use of the GUP 
eliminates the problem of an infinite temperature at the 
end of the evaporation process, which is clearly unphysi- 
cal, and leads directly to the prediction of the existence of 
black hole remnants ( 0, [H [HHH , HH). In references 
0, EH]: it has been shown that also the emission rate 
(erg/sec) is kept finite by the GUP mass-temperature 
formula, in contrast with an infinite output predicted by 
the Hawking formula. 



III. GOVERNING EQUATIONS 

In this section we will write down the basic equations 
which govern a system of black holes and radiation in the 
early universe. We will describe the evolution of a black 
hole mass as a balance of accretion and evaporation, as 
well as consider the dynamical behavior of a universe 
constituted by black holes and radiation. 
Then we shall derive a condition for a pre-inflation era 
dominated by matter, and the inflationary solutions for 
the equations of motion of the scale factor a(i), in both 
cases of pre-inflation matter, or radiation dominated eras. 



A. Emission rate equation 

In this subsection we will describe the evaporation be- 
havior of a micro black hole (in 4 dimensions) taking into 
account the GUP effects. 

In the present model we consider only photons or gravi- 
tons, nevertheless other kind of gauge or fermionic fields 
can be added in a straightforward way. 
Before writing down the emission rate equation, we re- 
view some delicate issues about greybody factors, emit- 
ted energy, and the Stefan-Boltzmann constant, in 4 di- 
mensions with the GUP. 

The presence of a GUP, i.e. of a minimal length, forces 
us to take into account the s queezing of the fundamental 
cell in momentum space (see [2ll . [26j-|29j ) . The squeezing 
results in a deformation of the usual Stefan-Boltzmann 
law. This deformation has to be considered, at least in 
principle, since we deal with micro black holes close to 
their final evaporation phase, where the predictions of 
the GUP are expected to differ noticeably from those of 
the Heisenberg principle. 

Due to the deformation of the Heisenberg fundamental 
inequality, 



hi U 2 
AxAp > - 1 + /3 Ap 2 



(16) 



the number of quantum states per momentum space vol- 
ume (or the invariant phase space volume) is 



dn x dn y dn z 



V 



dp x dp y dp z 



I Ai 2 



(17) 



Since p = hk, the number of quantum states (i.e. station- 
ary waves) in the volume V, with wave vector in [k, k-\-dk] 
is 



dN 



dn x dn y dn z 



V 



(27T) 



1 



(18) 



Since k — lo/c, the number of photons (or gravitons) with 
frequency within lo and lo + du in a volume V is given by 



dn 



V 

7 ~ ^ 



r» 



i 



h' 2 



3 e hu>/k B T _ l 



duj. (19) 



In case of a perfect black body (perfect emitter) we have 
for the greybody factor r 7 (w) = 1 for any lo. The depen- 
dence of T 7 (w) from the frequency lo is in general very 
complicated. It has been studied in many papers (for 4 
dimensional black holes see [30f . for emission of gravi- 
tons in 4 + n dimensions see [3l|), it is in some cases 
partially unknown, and in many cases can only be com- 
puted numerically. In the present model, we neglect the 
frequency dependence of T 7 , and therefore take the value 
T 7 := (r 7 (w)) averaged over all the frequencies. Thus, 
for the number of photons (or gravitons) in the interval 
(a;, lo + dto) in a volume V we write (in 4 dimensions) 



dn-, 



V 



[l + /3§(T i ) 2 ] 3 6 



huj/k B T _ I 



duo. (20) 



Obviously T 7 < 1 for a real non-ideal black body. 

The total energy of photons contained in a volume V (in 

4 dimensions) is then 



E^ OT {V) 



hio dn-. 



= r 7^§^r(4) C (4)A(/3,T), 



(21) 



where T(s) is the Euler Gamma function, £(s) is the Rie- 
mann Zeta function, and the function A((3, T) accounts 
for the cells' squeezing in momentum space, due to GUP. 
The function A(f3, T) can be formally written as 



A(/3,T) 



1 



P(4)C(4) Jo [1 + f3(2£ p k B Tx I he) 2 



■ dx 



and by this definition we have 

A(0, T) -> 1 for p 



0. 



Defining the Stefan-Boltzmann constant (in 4 dimen- 
sions) as 



CT3 



c T(4)C(4) 4 
3 ' ""' 



2 c 3 h 3 --B' 



the total energy can be written as 

, 3cr 3 



E^ OT (V) 



VT 4 A(/3,T). 



(23) 



(24) 



(22) 
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The energy dE radiated in photons (or gravitons) from 
the black hole, in a time dt, measured by the far observer, 
can be written as 



dE = r^V 3 T i A((3,T), 



(25) 



where V3 is the effective volume occupied by photons in 
the vicinity of the event horizon, 



V 3 = 47T Rg c dt. 



(26) 



Thus, fi nally, th e differential equation of the emission 
rate is [Hi Si SI 

-^- = l2nr i a 3 R 2 s T i A(f3,T). (27) 

where the minus sign indicates the loss of mass/energy. 
With the explicit definitions of 03, i?s, and using Planck 
variables m = M/M p = E/£ p , 6 = T/T p , r = t/t p 
(where £ p = ^ksT p and t p = £ p /c), we can rewrite the 
emission rate equation as 



dm 
~d7 



1 7r r 



7 m 2 Q A A(P,Q), 



15 



(28) 



where we used r(4)£(4) = 7r 4 /15 and 



tt 4 7o [1+4/36 2 ; 



dx 



(29) 



In the applications presented in the following sections, 
the GUP will be implemented by considering only the 
cutoff imposed on minimum masses and maximum tem- 
peratures. In other words, we mimic the cutoff effects 
of the GUP by simply stating that the micro black hole 
evaporation stops when T — T max or equivalently when 
M = M m j n , and otherwise using the "simpler" Hawking 
form of the mass-temperature relation. This is tanta- 
mount to choose A(j3,Q) ~ 1. We adopt this choice in 
order not to render the calculation too tedious, in par- 
ticular for those involving the nucleation rate of black 
holes and the emission rate in the presence of absorption 
terms. 



Absorption terms in the evolution equation for 
micro black holes 



In this section we consider the absorption of radiation 
by black holes. Therefore, we extend the emission rate 
equation (|28[1 in order to describe all the processes chang- 
ing the mass of a black hole. In principle, as our system 
consists of radiation and micro black holes, we should also 
take into account the absorption of micro black holes by 
other micro black holes. However, we will see later that 
the black hole density is low enough to neglect scatter- 
ing processes among black holes themselves, and thus we 



limit the analysis of absorption to the background radia- 
tion. The calculations are particularly inspired by Refs. 
[34U38j ]. Absorption terms will appear with a positive 
sign in Eq. (|28[1 . The general form of the absorption term 
will be 



dM 
~dt 



P 



(30) 



which contains the effective energy density p e ^ and the 
appropriate cross section a for the gravitational capture 
of relativistic particles in the background by the black 
hole. Since we want to consider relativistic background 
radiation, the effective energy density can be defined as 



P 



eff 



P + 3p(p) ■ 



(31) 



In the case of radiation with an equation of state param- 
eter w = |, this results in an effective energy density 
of 



P eff = Prad + 3^^- = 2p rad . 



(32) 



Thus, the absorption/accretion term for background 
radiation reads 



dM _ Orad 

'-'Prad ■ 

at C 



(33) 



Since the environment is supposed to be isotropic and 
homogeneous, the cross section for the absorption of rela- 
tivistic particles is proportional to the square of the black 
hole mass (35| . 



®rad 



'part 



27tt 



G 2 M 2 



(34) 



Note that a heuristic deduction of such cross section can 
be obtained directly from the spherical geometry of the 
black hole dM = 4ttR 2 s p rad dt/c = 16 tt G 2 M 2 p rad dt/c 5 . 

In Planckian units the equation for accretion terms 
reads 



dm 
~d7 



27nm 2 



n eff 



Ppl 



(35) 



where p p i := £ p /£ p is the Planck energy density. 

Then, a more complete differential equation for the 
evolution of the mass of a micro black hole can be given 
by 



dm 



15 



m 2 e 4 ^(/3,8) + 54 7rm 



2 Prad 

p P i 



(36) 



where we used the background equation of state through 
the effective energy density p e ^ = 2p rad . 
As stated before, for sake of simplicity we assume that 
the black hole evaporation evolves according to the stan- 
dard Hawking mass-temperature relation (I12[) . and thus 
we consider in Eq. (|36p the GUP correction function 
A(/3,Q) ~ 1. We shall keep in mind the cutoff on 
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mass/temperature predicted by the GUP, and put it in 
by hand whenever needed. 

Then, the differential equation for the evolution of mi- 
cro black hole mass/energy e can be written as 

de 1 o t r2 ^4 , 54ttG 2 £ 2 

-j£ = -127rr 7 (T3 R S T -\ p rad , (37) 

where e is the average energy content of a single black 
hole, e = Afc 2 . Using expression p3"|) for cr 3 , and i?s = 
2Ge/c 4 , we can write 



cfe C 2 



with 



3840 ttG 2 



54ttG 2 



(38) 



(39) 



C. Evolution equations 

1. Cosmological equation 

Given the standard RW metric (with Weinberg con- 
ventions but c^l) 



ds 2 



c 2 dt 2 + a 2 (t) 



1 - ki 



T dr 2 +r 2 dn 2 ) (40) 



where dil 2 = d9 2 + sin 2 Odtfi 2 , and the energy- momentum 
tensor of a perfect fluid 



Ty-v = (p + Pju^Uv + pg^ 



(41) 



where p is energy density and p is pressure, the (00) com- 
ponent of the Einstein equation reads 



ay kc 2 8ttG 
a 1 a 



1 3c 2 Pl 



(42) 



while from the (ii) components we have 



„fi fd\ 2 kc 2 8ttG , . 

2-+ - + — = -P- 43 

a \a J a z c z 

In our model, the energy density has contributions of 
radiation and matter, and can thus be written as p = 
Prad + Pmbh- For simplicity, and following Ref.Q, we 
now consider a fiat metric, i.e. k = 0. The equation is 
then written as 



87rG . . 

„ 2 (Prad + Pmbh)- 

O C 



(44) 



2. Evolution equations for p m bh and p ra d 

We suppose our system to consist of a "soup" of micro 
black holes and radiation. It is well known [see e.g. text- 
books by Weinberg or Landau] that for the description of 



the evolution of the energy densities p ra d and p m bh under 
the cosmic evolution of the scale factor a(t) usually the 
(0) component of the continuity equation is considered, 



G u 



(45) 



Here G° is a source-sink term that can appear especially 
in the description of reciprocally interacting subparts of 
the whole system, as we shall see in the next section. 
T^ v is the energy-momentum tensor of a perfect fluid. 
Specializing this equation to the RW metric, we obtain, 
for the global energy density p rad + p mbh , 



Prad + Pmbh + 4i? p ra d + 3i? Pmbh — G° 



(46) 



where H = a/a and G° is a possible source-sink term. We 
shall now compute accurately the form of the continuity 
equation for both the subsystems " radiation" and " black 
holes" , in particular the form of the source-sink term. 

The system we are investigating is a defined mixture of 
radiation and black holes, where the Hubble radius Rh 
of the universe contains a given fixed total number N of 
micro black holes, a given amount of radiation, and the 
only processes that can happen are exchanges of energy 
between the black holes and the surrounding radiation. 
As already mentioned before, in this phase no black hole 
merging, nor black hole nucleation, is supposed to happen. 

Let us first focus on the evolution equation for p m bh- 
Micro black holes are a particular type of dust: they can 
emit or absorb radiation. As a first step however, we 
suppose that micro black holes have a negligible interac- 
tion with radiation (i.e. we treat them as standard dust). 
Then the continuity equation, without any source term, 
can be written 



Pmbh + 3Hp m bh — \ 



H= a - 

a 



(47) 



This equation takes already into account the variations 
in mass/energy density due to the simple variation in 
volume. In fact, from (jT7|) we have 



and therefore 



Pmbh _ g / a 
Pmbh 



Pmbh 



(48) 



(49) 



where A is an integration constant. Hence we see that 
at any time t we should have p m bh(t)a(t) 3 = A, so the 
integration constant should be written as 



A = Pmbh(t c ) a(t c y 



(50) 



where t c is the time point when the constant A is deter- 
mined. It is a characteristic time for the onset of matter 
era, and will be investigated in section HVl Since a(t) is 
adimensional, A should have the dimensions of an energy 
density. 
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Conversely, considering black holes as dust grains of con- 
stant mass M, then the link between mass/energy density 
and volume can be immediately written as 



Pmbh 



Mc 2 N a{t c Y 



(51) 



where N is the total number of micro black hole in the 
volume a{t) 3 at any instant t > t c . As will be derived 
in section HVl N is considered to be constant, since no 
creation, merging, or complete evaporation of micro black 
holes are allowed after the time t c . All the quantities 
M, t Cl N, Rh^c), a(t c ) will be computed explicitly via 
numerical simulation in section IPV1 So we have 



. Mc 2 Na(t c ) 3 . _ 

Pmbh ~ 6 R H (t c ) 3 a 4 °~ 



Pmbh 



(52) 



which coincides with (|4"T)) . 

If now we suppose that also the mass/energy of the single 
black hole can change in time, then Eq. (|5Tj) reads 



Pmbh — 



e{t) N a(t c f 



R H {t c f a(t) s 



(53) 



with e(t) = M(t)c 2 , and this expression immediately sug- 
gests by derivation the correct source term in the conti- 
nuity equation: 



l',r,hh + 3 ( — ) p m bh — 



Na(t c ) 3 
R H (t c ) 3 a 3 



(54) 



Let us now focus on the equation for the radiation 
energy density p ra d- As first step, consider the variation 
of Prad due to presence of emitting/ absorbing micro black 
holes, when the system radiation/black holes is contained 
in a box of fixed volume. If de is the variation in a time 
dt of the energy content of a single black hole, and the 
box contains N black holes (all of the same mass), then 
the variation of the energy of the radiation in the box is 



dE = -Nde . 



(55) 



Since the volume of the box scales as V(t) — 
Rnitc) 3 a(t) 3 /a(t c ) 3 , then the variation of the radiation 
energy density is 



which means 



dprad — 



Prad — 



dE 



Nde 



V{t) V(t) 



Na{t c ) 3 



R H (t c 



(56) 



(57) 



This relation is true in the hypothesis of a fixed box. If in 
particular the black hole were inert (neither absorption 
nor emission) then e = and therefore p ra d — 0. 

In an expanding box, containing only radiation or just 
a few inert black holes (= dust grains), we know that 
from the continuity equation (|45[) we can write for p ra d 



Prad + ^Hprad = 



Prad 



B 
7^ 



(58) 



where B is an integration constant. 
Prom the previous two steps, it is then clear that, consid- 
ering both the expanding box and the emitting/ absorbing 
black holes, we can write globally for p ra d 



Prad 



Prad 



Na(t c ) 3 
R H {t c fa 3 



(59) 



We see that in the equations (|54|). (|59j) the term 
N a(t c ) 3 £/ Rnitc) 3 a? , which accounts for the emit- 
ting/absorbing activity by black holes, appears with op- 
posite signs, respectively. This is physically very plausi- 
ble since, if for example e > 0, then that term contributes 
to the accretion of black holes' masses, while exactly the 
same amount of energy is taken from the radiation sur- 
rounding the black holes. Coherently, we see that the 
global continuity equation for black holes and radiation 
combined reads 



Pmbh + 3 I - ) p mbh + prad + 4 ( - I p r ad = , (60) 



that is, it does not contain any source term. This is rea- 
sonable, since our system contains only black holes and 
radiation, and therefore the global energy content must 
be conserved (only diluted by the cosmic expansion rate 
H 7^ 0). Systems of equations where one term appears 
as a source in one equation, and as a sink in another, are 
quite common in physics and in cosmology. For example 
recent models dealing with the interaction between dark 
matter and dark energy display such features [39l ] . 



3. Complete set of equations 

We are now able to write down a set of equations that 
should describe, hopefully in a complete way, the pri- 
mordial "soup" of radiation and micro black holes, in a 
temporal interval ranging from the end of black hole pro- 
duction era (t — t c ) to the starting of inflation (t = ti n fi). 
Considering the equations (|55]). (l54")l . (p)9")) . pi)) , we can 
write the system (t is, as before, the comoving time) 



de 
lit 



C 



De Prad 



Pmbh + 3 ( — ) Pmbh 



Prad 

a 
a 



4 ( - ) Prad 



Na(t c ) 3 
R H (t c ) 3 a 3 
Na(t c ) 3 



R H (t c ) 3 a 3 



8irG . . 

„ 2 (Prad + Pmbh) 

o c 



(61) 



As we see, this is a system of 4 equations for the 4 un- 
knowns e(t), pmbh{t), Prad{t), a(t). This is a good sign for 
the closure and solvability of the system. However it is 
clear that this system is strongly coupled, and moreover 
nonlinear. Thus, to find an explicit solution is surely hard 
and perhaps impossible. Nevertheless, the system can 
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be studied in some physically meaningful situations (as 
for example when the micro black holes are very weakly 
interacting with radiation, with e ~ 0, when they essen- 
tially behave like dust). In these limits the equations can 
yield useful insights on the behavior of the scale factor 
a(t), which can be used (via a procedure similar to that of 
Ref. Q ) to compute the effects of this " soup" of radiation 
and micro black holes on the successive inflation era, and 
possibly on particular features of the power spectrum. 



D. Pre-inflation matter era 

We study here the regime just sketched at the end of 
the previous section, when micro black holes are very 
weekly interacting with radiation (e ~ 0). Our primor- 
dial "soup" is therefore composed by radiation and dust. 
In this approximation, the second and third equation of 
system (|61|) can be immediately integrated to give 



Pmbh o 



_ B 

Prad — J ■ 
Or 



(62) 



where the integration constants A and B have dimensions 
as energy densities, and can be written as 



.4 



Pmbh(t c ) d(t c ) 3 



B 



Prad(t r ) a{t r ) 4 . (63) 



Here t c and t r are the characteristic times for the onset of 
matter and radiation eras, respectively. They will both 
be explicitly specified in the next sections. Then the (00) 
equation of system (|6"Tj) reads 



8ttG 



A_ 

7^ 



(64) 



Equation (|64|) is separable, and can be integrated exactly. 
The solution obeying the boundary condition a(0) = for 
t = is 



(Aa{t) - 2B) ^/Aa{t) + B + 2BVb] = t (65) 



3^A 2 

where k = (87rG)/(3c 2 ). Using the binomial expansion 

vl/2 



(1 



, 1 1 2 

1+ r-r + 



we find, in the limit A — > or equivalently a(t) 
t->0, 



(66) 
for 

(67) 



which is the well know behavior of pure radiation era. 
In the other regime, when a or A are large, or equivalently 
B — » 0, wc find 



a 3/2 = t 

2 



(68) 



which is the the standard result for a matter dominated 
era. 



We can now wonder how much matter (micro black 
holes = dust) should be present in order to have a matter 
dominated universe before the beginning of the inflation. 
A condition for this can be easily derived by inspecting 
the exact solution (|65[) and considering its expansion as 




= t. 



(69) 



Clearly, the universe will be in a matter dominated era 
at the onset of inflation, namely for t = ti ~ 10 6 ip (the 
time when the temperature of the universe corresponds 
to the GUT energy scale), whenever the condition 



3 B 

2-a <<l 



(70) 



is satisfied. An even simpler derivation can be found by 
writing Eq. (|64p in the form 



A , 

= V 11 



B 
Aa 



(71) 



from which we read off that the evolution is matter domi- 
nated when B/(Aa) <§; 1. In section ITVl we shall compute 
explicitly via numerical simulation every step of the black 
hole nucleation phase, and the associated evolution of 
the pre-inflationary radiation and matter eras. We shall 
conclude that the above matter dominance condition is 
always met, even well before the onset of inflation. 



E. Inflationary solutions 

In our equations we now also take into account a con- 
stant vacuum energy, namely a cosmological constant. In 
this way, we shall be able to generate inflationary expo- 
nential solutions. Following Ref. J|, the (00) component 
of the Einstein equation in this case reads 



.4 

^3 



B 
7a 



C 



(72) 



The constant C in the Friedmann equation results from 
assuming a power law potential V{<j)) for the inflaton 
field. C and the potential are connected by 



v(4>) 



C-c x 



c 2 , 



(73) 



where C is the quasi de Sitter parameter in the Fried- 
mann equation, and c\ and ci are constants to be fixed 
from the inflationary model. In other words, C is mim- 
icking the potential for the inflaton field. The previous 
equation can be written as 



-V = K— (l — 

a) ~~ K a? V Aa A 



^1 («) 
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and under the matter era condition, B/(Aa) <C 1, it 
becomes 



>< 1 4- + c 



(75) 



Again, this equation is easily separable, and the solution 
obeying the boundary condition a(0) = for t = is 



1/3 



sinh I -V kC t 



2/3 



(76) 



which, for vanishing C, or small t, results in solution (1681) . 



a(t) 



2/3 



{nA) 1/3 £l\ 



(77) 



It is also useful to derive a condition for the onset of 
inflation. From Eq. (TT2")) we can obtain the sign of a 



a\ I A B 



and a > if 



_ A ( 2B\ A 



(82) 



(83) 



where we used the pre-inflation matter era condition 
B I (Aa) <C 1 . Therefore we shall be in a full inflationary 
era, a > 0, when 



> 



A 
2C 



1/3 



(84) 



while for large t^it exhibits an exponential (i.e. inflation- 
ary) behavior [401 ] . 



1/3 



exp 



(78) 



We can grasp an idea of the overall solution a(t) by nu- 
merically integrating Eg. (172")) . The evolution of the scale 
factor is shown in Fig [1] 



a(t) 




0.5 1 1.5 2 2.5 3 3.5 



FIG. 1: Diagram for a(t) versus t, in a model assuming sub- 
sequent radiation and matter eras before inflation. 



In case of radiation dominated pre-inflation era, i.e. no 
matter present (A = 0), equation (|T2"j) reads 



= /,■ ( £ + C) . (?<)! 



and the solution obeying the boundary condition a(0) 
for t = is 



a(t) = 



1/4 



sinh 



1/2 



which for small t or vanishing C is 

a{t) = V2 (kB) 1 / 4 t 1 ' 2 . 



(80) 



(81) 



IV. BLACK HOLE NUCLEATION: 
NUMERICAL SIMULATION 

In this section we numerically simulate the nucleation 
of micro black holes in pre-inflation era, and come up 
with a number of micro black hole remnants sufficient to 
make the universe pass from a radiation dominated to a 
matter dominated pre-inflation era. We will express ev- 
ery quantity in planckian units, meaning e.g. r = t/t p , 
= T/Tp, m = M/Mp. Thus every quantity is dimen- 
sionless. 

In 1982 Gross, Perry and Yaffe [ij investigated the stabil- 
ity of flat space and the arising gravitational instabilities, 
which might lead to singularities. They used the formal- 
ism of path integrals in a quantum version of Einstein's 
theory of gravity to analyze these gravitational fluctua- 
tions. As a concrete example, they took a box filled with 
thermal radiation to derive an expression for the prob- 
ability for the spontaneous formation of black holes out 
of the gravitational instabilities of spacetime. Two years 
later, Kapusta @ gave an alternative heuristic derivation 
of the nucleation rate, using the classical theory of nu- 
cleation during a thermodynamical phase transition. He 
reproduced the rate nearly completely with the classical 
approach considering the change in free energy of the sys- 
tem during the nucleation of a black hole, and completed 
the analogy by inserting by hand quantum corrections 
into his classically derived formula. See Appendix 1 for 
the explicit steps of such derivation. 
The nucleation rate for black holes reads fEq. (|209l) . Ap- 
pendix 1) 



rjv(e) 



8tt 



15 • 64tt 3 



■ e 



-167/45 e -l/167r6 2 



(85) 



where & is the temperature of the universe (the thermal 
bath), expressed in Planck units, and at the same time 
the temperature of the nucleated black holes, connected 
to their mass m by 



e 



i 



47rm 



(86) 
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At a given temperature 9, all the black holes created 
according to this nucleation rate will have mass m. 

As stated before, the pre-inflation era is supposed to 
take place from the Planck time t p ~ 10~ 43 s to the onset 
of inflation, ti n f ~ 10~ 37 s, when the temperature of the 
universe has reached the GUT energy scale. At very early 
times, when > 0* ~ 1/(4%), the nucleation probabil- 
ity is very high, but does not lead to black hole formation 
as it is forbidden by the GUP to create smaller than the 
Planck mass black holes (m m i n — VC M p where £ ~ 1). 
For temperatures above 0* ~ -h- production of small 
black holes is not possible. So at least for this very early 
time, the universe is simply a chaotic hot sphere that 
we suppose to be filled with primordial radiation, follow- 
ing the approach of Refs.[2|, |3j. There might be regions 
with larger density than others, but the conditions are 
too chaotic to allow formation of stable objects like black 
holes. The universe can thus be assumed to be radiation- 
dominated at the beginning, and will migrate to being 
matter-dominated at a later time, when black holes are 
starting to be formed. Considering an adiabatically ex- 
panding universe, we can write T(t)o(t) = T p a(t p ), and 
since during the radiation era the scale factor evolves like 
a r (r) ~ a(t J ,)r 1 / 2 , and we choose a(t p ) = 1, we have 



during the time when the universe is dominated by radia- 
tion. For a matter-dominated universe, the temperature 
evolution is analogously given by 




A. State of the universe 

The parameters in the Friedmann equation contain- 
ing matter and radiation energy densities are defined as 
(Eq.©) 

A = p m (Tc) ■ a 3 (r c ), (89a) 
B = Pr {T p )-a 4 {T p ), (89b) 

and their dimension has to be energy per volume, as the 
scale factor is a dimensionless quantity. 
At the Planck time, the universe is in a radiation- 
dominated stage, and it is reasonable to assume that 
Pr{r p ) = p p . Assuming a(r p ) = 1 (which is an unconven- 
tional, but convenient choice, and will be converted to 
the common notion of a(T to day) = 1 later), the radiation 
parameter can be fixed as B = 1, expressed in Planck 
units. For the matter parameter, we have to choose a 
time r c when black holes are starting to be nucleated, 
and calculate p m (T c ). This time and the parameter A 
will be derived in the next subsections via numerical sim- 
ulations. 



To calculate the time of the transition from radiation- to 
matter-dominance (which evidently takes place after r c ), 
we consider the ratio 

R(r) = = ^iM . r 2. (90 ) 

Pr(T) Pp 

At the very beginning there is only radiation. So for a 
while p m { T ) — 0, the evolution of P t (t) is driven by radi- 
ation, o(r) ~ t 1 / 2 , and therefore p r ~ P p /t 2 . This is cor- 
rect at least before black holes are created, whereas after 
the onset of nucleation, a(r) evolves in a more complex 
manner dictated by equations (|64l l65|) . As nucleation 
starts, p m { T ) grows, and some time later R(t) crosses 1. 

R(t) can only be used for qualitative statements at the 
beginning of nucleation, since, in the way it is defined, 
it is rigorously valid only until the start of matter nu- 
cleation, and it does not contain any information about 
the time-development of the scale factor according to the 
full Friedmann equation, when matter and radiation are 
both present. So it should only be applied during short 
time spans, just after the nucleation starts, when the 
scale factor doesn't change significantly. The condition 
3B/(2Aa) <C 1 is the only significant criterion to fully 
determine the radiation to matter transition of the uni- 
verse at later times, when matter and radiation are both 
present. 

B. Nucleation Process 

Using Eq. (|87|) , in radiation era we can write the nu- 
cleation rate as a function of time as 

Tn At) = — ^— . T 167/90 e -r/16 W _ (91) 

15 ■ oir z 

The temperature of the Universe is linked not only to 
time, but also to the mass of the nucleated black hole, 
since a black hole, at the instant of its creation, is in 
thermal equilibrium with the rest of the Universe. In 
fact 

= -4— and r = —4 => T r = 16ir 2 m 2 . (92) 

It should be explicitly noted that the relation (|9"2"|) does 
not express the time evolution of the mass of one black 
hole, but the dependence of the initial mass of the nu- 
cleated black holes on time. The evolution with time of 
the black hole mass, due to evaporation/accretion pro- 
cesses, is given by relation (|36p . while equation ex- 
presses the evolution with time of the masses of the black 
holes at the instant of their creation. Therefore the cutoff 
from the Generalized Uncertainty Principle, which gives 
a minimum mass m ~ 1, can be translated also in terms 
of time as t* = 167r 2 ~ 158. This can be seen in Fig. [2]- 
the curve is truncated at r + (vertical line), which corre- 
sponds to the cutoff at about one Planck mass. This re- 
lation also implies that the black holes nucleated at later 
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FIG. 2: Nucleation rate, Tn,t{t) over time r. 



times have larger masses, according to Eq. (|92)) . while the 
probability of their formation decreases. 

From the onset of nucleation the number of black holes 
nucleated each Planck time per Planck volume is given by 
the rate (f9Tj) , but we have to monitor closely the over- 
all state of the universe. When the ratio R(t) crosses 
unity, then the universe changes to a matter-dominated 
stage, in which case the nucleation rate is no longer given 
by Eq. (f9"Tj) . As soon as the phase transition happens, the 
nucleation rate has to be given in terms of the tempera- 
ture in a matter-dominated universe (0 m = t~ 2//3 ): 

T N , m (r) = • r 334 / 135 • e- ^' ^. (93) 

In order to get the number of black holes nucleated with 
time, we have to do a time integration over the nucleation 
rate, multiplied with the volume of the universe in each 
time point. Considering that the Hubble radius is defined 
as Rh — c/H, and that for power law expansion rates 
(a(t) ~ t a ) we have Rh ~ ct, then, taking as initial 
condition Rh(t p ) — £ p , we can write (in units of Planck 
length) Rh(t) = r. Therefore (since l z v — V p = 1) 

4-7T 

Vh(t) = y r 3 . (94) 

But before we go on to calculate the black hole number, 
we take into account another principle, which will give 
a more stringent cutoff on the nucleation of black holes 
than the GUP does. 



C. A Cutoff from the Holographic Principle 

The holographic principle [4l[ places a limit on the 
information content, or entropy content, in a certain re- 
gion of space-time. Quantitatively it states that (in units 
where ks = £ p = 1) 

S[L(B)} < ^1, (95) 



where L(B) is a so-called light sheet, which defines a cer- 
tain region of space-time B, and A(B) is the codimension 
1 boundary of that region. For our situation, applying 
the holographic principle simply means that the total en- 
tropy contained in the Hubble sphere cannot exceed the 
entropy of a black hole of size equal to the Hubble sphere, 
which represents the maximum entropy that can be held 
in that spacetime region, as black holes are the most en- 
tropic objects. According to Refs. [Ulllij], the expression 
for the entropy of a black hole is 

Abh 

^bh — — |-j 

and so the entropy of a black hole of the size of the Hubble 
sphere (HS) is given by 

S hs (t) = ^=7tR h (t). (96) 

This can be used to define a cutoff for the nucleation rate. 
We demand that at no time point in the evolution of the 
universe the entropy of the black holes can exceed the 
total entropy that Hubble sphere can maximally hold, 
and then we try to find a time point r c , from which this 
condition is fulfilled. If this condition is violated in the 
course of black hole production, then it is simply not 
allowed to create black holes. The entropy of the black 
holes within the Hubble sphere is 

S bh (r) = jy N Ar') ■ nr 2 s (r') ■ dr' , (97) 

where the Schwarzschild radius of the black holes is r s = 
m = t 1 / 2 /(47t), from Eq. ([9^|) . and we require that at all 
times 

Sbh(r) < S H s(t)- 

Again, since the universe has started off in a radiation- 
dominated era, and the black hole nucleation has not yet 
begun, the rate as a function of time is defined using 
r = t -1 / 2 . It is yet unknown when the production of 
black holes will change the state of the universe to be 
matter-dominated, but to find out the time of nucleation 
start, we are for now bound to the assumption of radi- 
ation dominance and thus will use Tpf ir (r) in the above 
expression for the entropy. 

The cutoff time turns out to be r c = 993, as can be seen 
in Fig-El which shows the entropy of the Hubble sphere 
(dashed line) and that of the black holes (full line) as a 
function of time. 

With t c being determined, we can calculate the number 
of black holes according to 

N bh (T) = ^ r N , r (T) ■ (yT 3 )dT', (98) 

and the matter density as 

Pm(r) = 4 ^ rjv,r(r') ■ m(r') • r' 3 dr'. (99) 
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FIG. 3: The cutoff as defined by the holographic principle at 
t = 993: the curves touch, but the entropy of the black holes 
(full line) never exceeds the maximum entropy that can be 
held by the Hubble sphere (dashed line). 



FIG. 4: Ratio R(t) of matter to radiation density in the early 
universe, after the nucleation starts at r = 993. 



However, in order to correctly calculate the black hole 
number and the mass density, we need to know when the 
universe will migrate from being radiation-dominated to 
being matter-dominated. We know that black hole nucle- 
ation starts when the universe is in a radiation-dominated 
state. For the very beginning of black hole nucleation, 
for a time span of at least one Planck time, we have 
to use the rate rjv r to produce black holes. After one 
Planck time, we have to check whether the universe is 
still in radiation- or already in matter-dominated stage, 
and take the according production and expansion rates to 
follow the black hole production. A new parameter d r is 
introduced, which denotes the duration of this radiation- 
dominated phase of black hole nucleation, until the uni- 
verse reaches a matter-dominated state. The matter den- 
sity of the black holes can be expressed like 

1 l-Tc+dr 

Pm(r) = ^3 / rjv,r(r') • m r (<r') • T a dT 



1 



IV, m (T') • m ro (r') • r'W, (100) 



T c +d r 



where m r = t 1 / 2 /(47t) and m m — t 2 ^ 3 /(47t). We have 
separated integrals for the two separate stages of black 
hole production. We can now vary the parameter d r and 
investigate the ratio we are looking for, (here p v = 1) 



R = p m (r) ■ t 2 



(101) 



to determine at which point the dominant substance in 
the universe changes. 

As can be seen from Fig.[H it turns out that the universe 
very quickly reaches a matter-dominated stage after the 
onset of black hole production. After only a few Planck 
times, the ratio is clearly above unity, and so we can 
safely assume that the universe is in a matter-dominated 
stage at about 3 — 4 Planck times after the black hole 
nucleation has begun. 



With this assumption the estimated number of primor- 
dial black holes is given by 

A(r m/ ) ~ 10 4 . 

The masses of the black holes nucleated between 
t c = 993 and r = 996, according to the mass-time- 
relation Eq. (|9"21 , are 

m(r c ) ~ 2.5. 

We can now also calculate and plot the mass den- 
sity, Eq. (|100[) in evolution with time. It is shown 
in Fig.E 
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FIG. 5: The mass density of the universe as a function of 
time. 



D. Collision Rate and Black Hole Thermodynamics 

What is left to do to justify the assumption of a matter- 
dominated era due to the existence of black holes is to 
analyze the results in the context of black hole collisions 
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and merging. The collision rate of black holes can be 
given by starting from a general definition of a scattering 
rate. We consider black holes of mass m(r c ), and their 
velocity to be determined by Brownian motion, therefore 

v = \J jj. Being n the number of black holes per unit 
volume, and a the scattering cross section, we arrive at 



^coii = nav 

= n bh ■ 47rrf • c\l — 
V m 

N bh (r) 
V h (t) 



1/3 



(102) 



where in the last line we used r s = m and 6 m = r~ 2 / 3 , 
with c = t p = 1. The rate can be seen in Fig. [51 together 
with a dashed line denoting the Hubble expansion rate. 
The collision processes are effective as long as the rate 
is higher than the Hubble expansion rate - otherwise the 
reaction decouples, and collisions stop because the ex- 
pansion of the universe is faster than the time between 
collisions. As can be seen clearly from the plot, the colli- 
sion rate is at all times lower than the Hubble rate, and 
thus no collisions or mergings of black holes take place. 



r.„u(rt,H(r)[io- 4 ] 
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FIG. 6: The collision rate of black holes (full line) in develop- 
ment with time, together with the Hubble rate (dashed line). 



Since a(t p ) — 1, we can write a(r) ~ 10 1 / 2 t 2 / 3 and 
B p r {t p )a*{t p ) 100 Pp 



Prad(T) - a4(r) - q4(t) - r8/3 

So the differential equation (f36|) becomes 

dm 1 5400 7r m 2 

~dr~ " ~480 7rm 2 + ^ 



(104) 



(105) 



which can be numerically integrated with the initial con- 
dition 



m(r = 998) ~ 2.5 



(106) 



The numerical integration confirms that, despite the ab- 
sorption term, the black hole evaporates completely in 
a time T(m) ~ 17086 t. p . Actually, without the ab- 
sorption term, the lifetime would be about 50% less, 
just T(m) = 7854 t p . In any case, by a cosmic time 
of t ~ 2 • 10 4 t p , well before the inflation starts, all the 
black holes have evaporated down to the Planck size rem- 
nants predicted by the GUP. In Fig. [7] the two lifetime 
diagrams, one for evaporating/absorbing black holes, the 
other for evaporating only black holes, are compared. 
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FIG. 7: Black holes lifetime diagrams: the upper (red) for 
an evaporating/absorbing black hole; the lower (green) for an 
evaporating only black hole. 



We can also do some estimations of the black hole ther- 
modynamics. As we have seen that the black holes are 
not colliding with each other and don't have a chance 
to merge, we know that the black hole masses can only 
change by accretion or evaporation. This process is de- 
scribed by the differential equation (|36l) . In our spe- 
cific case we consider, as usual, the correction function 
A(f3, 0) ~ 1 and the greybody factor T 7 = 1. 
From t c ~ 1000 to the onset of inflation the Universe is 
matter dominated, therefore we can write for the scale 
factor 



a(r) 



a(tp)r 
a(t p ) t c 



1/2 



1 < T < T r ~ 10 3 t K 



"1/6 _2/3 



(103) 



T r <r< T, 



infl 



10 b i„ 



E. Fixing the Friedmann Equation 

In this subsection we come back to the original goal of 
this part: to determine the numerical size of the param- 
eters in the Friedmann equation that correctly describe 
the universe developed in our model. We have already 
settled B — 1 by simply assuming that the density of ra- 
diation at the Planck time was equal to the Planck den- 
sity. For the matter component, we know that about 10 
black holes are nucleated during a short period around 
t ~ 10 3 in the pre-inflation era. A is given by Eq. (|89al) . 
The matter density can be simply estimated by 



Pm{T c ) 



10 4 black holes 10 4 v / Ce P 



10 9 K 
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where m min = ^/C 1S the minimum mass predicted by 
the GUP. The scale factor is chosen as a(t p ) — 1, and 
then evolves in radiation-dominance until r c ~ 10 3 ; thus 
a(r c ) ~ Tc^ 2 ~ 10 3 / 2 . Putting these together, the final 
result for the matter component is 



.4 



io 4 ^ 

10 9 K 



10 



9/2 



10 



-1/2 



(107) 



As stated before (see Sections III A[ III Cp . we consider 
£ ~ 1. With this number, we can confirm the actual 
validity of the matter-dominance condition for the whole 
era r c < r < Tj n //. In fact in this era, according to 
Eq. (|103p . the scale factor evolves over the values 10 3 / 2 < 
a(r) < 10 7 / 2 , and therefore the matter-dominance ratio 
spans the values 



B 



2 A -air) ~ 10" 1 / 2 • o(t) 



10" 



10~ 3 < 1 



which confirms that the Universe is in matter era from 
t c to the onset of inflation. 

With all this information, we are fully equipped to com- 
mence the calculations of the primordial power spectrum. 
It should be noted that C is not yet determined - this will 
be taken care of in the next sections. 



F. The Scenario without the GUP 

What if the nucleated black holes wouldn't be pro- 
tected by the GUP from evaporating completely but 
would just transfer their energy to radiation filling the 
universe? There would be no era of matter-dominance 
before inflation, but there would only be radiation 
until the onset of inflation. Or would it? In the 
previous chapter, we have found out that the black 
hole nucleation is very powerful and the universe is put 
into a matter-dominated stage nearly immediately. If 
the GUP is assumed, the evaporation stops when the 
mass of the black hole reaches Planck mass, and the 
Universe stands in matter era until the onset of the 
inflation. On the contrary, if the GUP is not valid, 
the black holes, once created, can evaporate down to 
zero mass. The universe remains in a matter-dominated 
phase for the time of evaporation, but when all black 
holes have vanished, radiation is the dominant species 
again, until inflation starts. That means, there are three 
stages during the pre-inflationary era, and the evolution 
and growth behavior of the scale factor is changing at 
three transition points in time. The parameters of the 
Friedmann equation have to be evaluated anew, and 
this leads to a new differential equation for the field 
perturbation and a new fc(a)-relation (see Sec IV A [I . 
We can summarize the time frame of this scenario as 
follows. The universe starts out radiation-dominated 
and remains so until the start of black hole nucleation 
at time r c ~ 10 3 , where the scale factor has expanded 
to about a(r c ) = 10 3 / 2 . The subsequent matter phase is 



short and only lasts from r c to r r ~ 2 • 10 4 , as we know 
from the previous calculations that the nucleated black 
holes have an approximate life span of 2 • 10 4 i p . During 
this time span the scale factor expands for a factor of 
10 2 / 3 . Again follows a period of radiation dominance, 
which lasts from r r to Ti n t. The temperature at the 
onset of inflation here again has to match the GUT 
energy scale, and for this reason inflation is starting a 
little later, at r, n / ~ 10 s . From r r to Tj n / there are 
four orders of magnitude in time, and so the scale factor 
expands for a factor of 10 2 . 

This is the situation we are facing when the GUP is 
not included in the theory. However, there is a third 
possibility that can be considered. 



G. The Scenario without any black hole 

After the previous subsection, it is also justified to ask 
for the case when there are no black holes at all, and no 
black hole nucleation. What happens if we basically ne- 
glect all the black holes nucleation processes in the early 
universe, and simply assume that the pre-inflation era 
was completely radiation-dominated? The results from 
this scenario will then actually correspond to the inves- 
tigations done before (e.g. @, yj]), completely cutting out 
the possibility of black hole formation. In those previous 
works, a suppression of the CMB power spectrum for low 
multipole modes, resulting from a suppression in the pri- 
mordial power spectrum, was found. 
Under these assumptions, the universe is in a radiation- 
dominated stage from the Planck time until the temper- 
ature reached GUT energy scales, at about r = 10 7 i p . 
During this time span, the scale factor can expand from 
one Planck length to 



10 7 / 2 



(108) 



The power spectra resulting from the scenario introduced 
here are computed in section IV El 



V. SCALAR FIELD FLUCTUATIONS ON AN 
EVOLVING BACKGROUND - THE 
PRIMORDIAL POWER SPECTRUM 

In this section we study the quantum fluctuations of a 
field living in a universe whose background evolves with 
a given scale factor a(t). To facilitate comparison with 
the case of a pre-inflation radiation dominated universe 
already considered in Ref.Q, we choose the same con- 
vention on the metric, namely 



ds 2 = g^dy^dy" = dt 2 - a(t) 2 dy 2 



(109) 



Here we set c = 1, and we choose a flat metric (for 
simplicity, and since we deal with an almost flat uni- 
verse). We consider a zero mass scalar field $(i, y) and we 



15 



perturb the field around its classical expectation value, 
Q(t, y) = $ (t, y) + (p(t,y). The equation of motion for 
the scalar field perturbation then reads 



where 



Uip{t, y) = 0, 



(110) 



Another possible strategy to solve Eq. (1113p is to change 
the independent variable from t to a. Then 

d da d . d ■ . ,. , „ s 

dt = TtTa= a da-> = ^ 

and the equation for <pk (a) reads 



□ 



(in) 



In the applications, the a(t) of Eq. ([76|) . or of Eq. (|80l) . will 

be used in the metric g^ v . 

Taking the Fourier transform for ip 



ip{t,y) = I cj> k (t) e-' lk -y dk, (112) 
we obtain from (|110[) the equation of motion for <fik(t), 



Mt) + 3-Mt) + (% ) fafjt) = 0. (1 I.:!) 



In order to solve Eq. ([113[) . different strategics are cus- 
tomarily used. 

One possibility is to introduce the conformal time r\ de- 
fined by the relation 



dt 

dq = -7a > 
a(t) 



which implies 



drj d 



1 d 



dt dt drj a(rj) dr\ 
Then the equation for 4>k{v) becomes 



<t>Uv) + 2-<P'k(v) + k 2 Mv) = 0, 
a 



(114) 



(115) 



(116) 



where a "prime" indicates a derivative with respect to rj. 
Applying now a well known general procedure, valid for 
any second order differential equation, we can make the 
first derivative disappear. In fact, defining 



with 



we find 



P(v) = exp 



p(v) 



a(r))' 



and henceforth the equation for Vk{rj) reads 



4 



k 2 v k = . 

a 



(117) 



(118) 



(119) 



(120) 



• 2i« 

1 





d 2 \ 






(- 


a J 


<& + 





b k = 0, (122) 



where of course d and d have to be expressed as functions 
of a. In the following, we shall use indifferently both 
procedures. 



A. The re-entering k-modes 

The k-modes which left the horizon at or shortly af- 
ter the onset of inflation are those that are just now 
re-entering our Hubble radius, and they represent the 
largest modes of fluctuations currently observable, having 
a size comparable with that of the visible universe. When 
we solve equation (|113[) . or equivalently Eqs. (|120l) . 
(I122p . we fix a particular value of the parameter k, i.e. 
a particular mode, and we obtain the solution for this 
mode, introducing suitable boundary conditions. We are 
concerned with modes that leave the horizon just at or 
shortly after the beginning of inflation, because they cor- 
respond to the largest scales observable today, and they 
could bring imprints of a possible pre inflation era. A prc- 
inflationary era affecting these modes could thus explain 
the anomaly of the quadrupole moment of the CMB to- 
day. Simple geometric considerations will help us to find 
a relation between k and the scale factor. By compar- 
ing the FRW metric (|109|) . which is written in comoving 
coordinates, with the standard euclidean physical metric 
dt 2 — dY 2 hys , we infer the relation between comoving and 
physical coordinates 



dY phys = a(t)dy c 



(123) 



Recalling the Hubble law, v — Hd, we get the physical 
Hubble radius as c = HRh p , and therefore the comoving 
Hubble radius as 



R 



H,c 



Rh,p 



C 

aH' 



(124) 



Now, the wavelength of a perturbation crossing the hori- 
zon at a time t c (the largest visible perturbation) should 
be given by the relation 



which means, in comoving coordinates, 

Ar, 



A r 



a 



a c H c 



(125) 



(126) 
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Therefore a particular comoving fc-mode at the time of 
its horizon crossing obeys the relation 



(127) 



Since we chose units where c = 1, and moreover tt/2 ~ 1, 
we can finally write 



aH 



(128) 



as the horizon crossing condition. Besides, reminding 
that H = a/a, we have for a pre-inflation radiation era 



k = aJ~n I —r + C 



1/2 



while, for a pre-inflation matter era, 



k = aJTi [ + C 



1/2 



(129) 



(130) 



In FigJSl we give a plot of the (k, a) relation in the case 
of pre-inflation matter era, with constants chosen arbi- 
trarily as A = C = 1. 




FIG. 8: Diagram for k versus a (full line), in pre-inflation 
matter era 



B. Scalar Field equations 

In this subsection we specialize Eq. (|122j) to particular 
a(t) solutions. For the pre inflation matter era case, using 
Eq. (I75|) or (|76p to compute a, a, we obtain 



1 /4Ca 3 



7-4 



Ca 3 + A 



na(Ca 3 + A) 



<Pk 



0. 

(131) 



In the case of pre inflation radiation era we have, using 
Eq. flU, 



The general plan of our work is to solve equations (|131l) , 
(|132l) . regarding k as a parameter, to obtain <j>(a, k). Al- 
ternatively, we can solve Eq. (|120[) and get (j>(r},k), and 
then, since a = a(rj) and 77 = 77(a), make a substitu- 
tion to 4>(a, k). Once the solution <p(a, k) is available, we 
can compute the power spectrum P(k) of the quantum 
fluctuations of the field $. Since the power spectrum is 
usually given as a function of k only, will be necessary to 
express <j> as a function of k only, and this can be done 
through the relations (|129l) . (I130p . The final scope is to 
obtain the function 



P{k) = fc 3 |<?(>(a(fc),fc)| 5 



(133) 



which represents the primordial power spectrum of the 
fluctuations (perturbations) of the field $. With P(k) 
then we feed - as in our case - the CMBFAST code to 
generate the CMB anisotropy power spectra. 

Unfortunately, analytical solutions of the full equa- 
tions (IT5T|) . (|T3"2"|) . (|T2"0")l . when we deal with functions a(t) 
given in (|76|) . (|80)l . are not easily expressible in closed 
form. A power series solution is in general at hand, but 
still not much more comfortable than the mere numerical 
one. However, luckily enough, it turns out that analyti- 
cal solutions of the second order approximations of equa- 
tions (|131l) . (|132l) can be obtained by means of the WKB 
method, and these solutions are strongly corroborated by 
numerical insights. 

Of course, when solving a differential equation, we need 
boundary conditions in order to fix the solution explic- 
itly. From the mathematical point of view, a boundary 
condition can be put anywhere in the realm of definition 
of the solution. From the physical point of view, it is wise 
to put boundary conditions in regions where the physics 
is reasonably well known. Since we are almost sure that 
inflation happened, while we know little about a possible 
pre-inflation era, we choose to put our boundary condi- 
tions in the full inflationary era. This means that, for 
large a, or equivalently, for large k (see equations (I129p . 
([130])), well after the beginning of inflation, the field 4> 
must generate an almost scale invariant, i.e. flat, primor- 
dial power spectrum P{k). To be more precise, the latest 
observations have shown that the primordial power spec- 
trum is not exactly flat, but slig htly tilted. The newest 
analysis of the WMAP data [45| indicate a form like 



P(k) 



with 



n s = 0.963 ± 0.012 (68% CL) . 
Therefore the field 6 must behave as 



\ct>(a(k),k)\ 



F/2 



(134) 



(135) 



(136) 



Ca A 



Ca A + B 



1 



n{Ca 4 + B) 



= 0. f° r large k. 

This condition will allow us to fix properly the arbitrary 
(132) constants in our solutions. 
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C. Analytical investigations: First order 
approximation 

To begin with, let's figure out the solution of our dif- 
ferential equation for very large a, or k. We use, as first 
example, the equation for Vk(rj), (|120|) . We need the ex- 
pression of the conformal time 77 at large a, i.e. large 
cosmic time t. In fact, for t — > 00, we have from Eq. (j?6")) 



a m (t) 



C 



1/3 



i_ \ exp (y^t) , (137) 



sinh I - v kC t 



1 2/3 



and for the radiation solution, Eq. (|80| . 



a r (t) 



D 

C 

B 

AC 



V4 r 



1/4 



sinh ( 2VkC t 



1/2 



exp 



(138) 



Note that for large t the solutions are essentially the 
same, as it should be, since inflation "washes out" ev- 
erything (actually, almost everything, as we shall see). 
From dr\ = dt/a(t), we can express a(t) with the confor- 
mal time 77, 



0(77) 



1 1 



1 



\v\ 



'kC t) 



(139) 



with t) < 0, and this result holds, in the limit t — > 00, for 
both solutions a r (t) or a m (t). Thus, the equation for Vk 
reads 



Ufe = 0. 



(140) 



The general solution in terms of Bessel functions is 



Mv) = VW\ [ci(fc) J 3 / 2 (k\7)\) + c 2 (k) J_ 3/2 (fe|77|)] 



(141) 



where 



J-3/2(x) 



2 1 
2 1 



Sill X 



cos a; 



/ 1 / cos 2; \ 
= ( sin x 

7T Jx \ X ) 



with 



Then 



x = k\q\ 



' kC a 



(-2 



2 1 

(*) ( 



Cl(fc) 

cos a; 



cos a; 



X 

smi 



(142) 



Let's have a look at the behavior of the solution for a 
00, which means \rj\ — ¥ 0, x — >■ 0. Then 



ci(fc) x 2 + c 2 (fc) [I + £-?L 



(143) 



Since x ^ 
4>k{a) 
Ufe(a) 



we have, for any given fc, 



*(^+<*(*)(l + £i 



(144) 



1 

8aA 



If we require, as usual, 4>k(a) — » for a — > 00, then this 
necessarily implies 02(h) = . So finally the solution is 



«k(»?) = [ci(fc)J 3 /a(fe|»?|)] 



and 



4>k(a) 



Vk(a) 



(145) 
(146) 



IkJC 



■ sin 



'nCa 



cos 

a 



' nCa 



This is the solution which describes the field </>(a, fc) in 
full inflationary era. As we have seen, this solution holds 
for both cases of pre-inflation matter and pre-inflation 
radiation era. The arbitrary integration constant ci(fc) 
still has to be fixed, and this can be done by enforc- 
ing the boundary condition in inflationary era, namely 
Eq. (|136j ). Noticing that, for large a, from both equa- 
tions (|T2"9" J) .([T3"0" ]) . we have 



'kC, 



(147) 



and therefore 

0(o(A),fc) 



^ [sin(l)-cos(l)] (148) 



then this means, from (|136|) . 

d(jfe) = k^ n ' 



(149) 



The solution (|146l) will be employed as boundary condi- 
tion in inflationary era for the numerical integration of 
equations fl3"I|) . (fl%2"|) . 

Finally, we note that we can start from equations (|131l) , 
(|132[) for 4>k (a) , instead of equation (|120j) for Vk (77) . Tak- 
ing the first order approximation for a — > 00 in the coef- 
ficients, we get, from both equations (|131[) . (|132[) 



*k = 0. 



This has the solution 



<Ma,fc) = ^Uc 2 (fc), 



(150) 



(151) 



which, for large a, when fc ~ a V kC, can be matched 
with the "almost flat spectrum" boundary condition 

CDSJ) if 



d(fc) ~ fc I + ^ n =- 1 ) 



c 2 (fc)=0. (152) 
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D. Analytical investigations: Second order 
approximation 

Obviously, a solution like (I146[) cannot be trusted when 
we look for information about the behavior of P(k) at 
low fc-modes, since by construction it is built by taking 
into account the first order approximation only, and re- 
quiring that it reproduces the almost flat spectrum for 
high k's. To have reliable analytical insights on the low 
k behavior of P{k), one should go to the second order 
approximation. This can be done in principle by starting 
from equation fjl 20[) for Vk, but it would involve a quite 
complicated construction of the conformal time rj, and of 
the function a{rj), in terms of Jacobian elliptic functions 
(see Appendix 2). A more straightforward path turns 
out to be to start directly from equations (|131[) . (|132[) 
and develop the coefficients to the second order in a for 
a — > oo. 

For the pre- inflation matter era case, Eq. (jl3ip . we have 



4 3A 1 _ / 1 
+ O -= 



a 



2C a 4 



1 



k Ca 4 



0[- 

a' 



b k = 0. (153) 



For the pre-inflation radiation era case, Eq. (|132p . we 
have 



4 2B 1 
a ~ ~~C~^ 

k 2 1 



<> ' 4 ) ) 



k Ca 4 



O 



h = 0. (154) 



Following the usual well know method to get rid of the 
first derivative, we write 



(a) = v(a)p(a), 



with 



p(a) = exp 



-- / f(a)da 



(155) 



(156) 



where f(a) is the coefficient of cj>'. Then the equation for 
v(a) reads 



v"(a) + [ t - + f--+g)v(a) = 0, (157) 



P 



P 



where g(a) is the coefficient of 0. 

To an equation of the form (| 157[) we can apply the WKB 
method. In what follows, we specify the main steps of 
the argument for the matter era case, and report only 
the final result for the radiation era case. 
For Eq. (|153p . pre-inflation matter era, we find 



1 / A 1 

pw = — ex P -77;^ 

a z \ 4C a 6 



(158) 



and thus 



Setting 



2 
~~ 2 



9 A 2 1 . 

- c — 2 -e)v = ® (159) 



F(a) = i 



9 A 2 1 



kCci 4 16 C 2 a 8 



1/2 



(160) 



c + (fc)exp i / F(a)da 



then the WKB ansatz suggests as solution for v(a) 
v(a) — . ^ 

+ c_(fc)exp^— i J F(a)iia 

Writing 



(161) 



G{a) 



F{a)da 



(162) 



we see that, for example, for a — > 00 

G{a) ~ iv^loga. (163) 

So, the WKB solution of Ea. (|153p explicitly reads 

(j>k{a) = Vk(a)p(a) = (164) 
2\^C 2 [c + {k)e lG ^ +c_(fc)e- 4G ( a )] e"™/ 4 
[32nC 2 a e - 16k 2 Ca 4 + 9kA 2 } 1/4 ■ exp[A/ {4Ca 3 )} ' 

We can always find agreement with the boundary condi- 
tion (|T36"1) . namely \<f>\ - fc-3/2+(n 3 -i)/2 for large a Qr k ^ 

by defining the arbitrary constants c±(k) accordingly. In 
fact for large a, a ~ fc/V kC, and we can choose c±(k) in 
a way that 



c + (fc)e 



»G?(o(fc)) 



c_(fc)e" 



-iG(a(fc)) 



fc(n.-l)/a ( 165 ) 



so that the condition (|136|) is fulfilled. The expression of 
4>(a(k),k) when a(fc) ~ kj\J nC (first approximation) is 



0(a(fc),fc) 



2 v^C2 fc K-i)/2 e -»/4 



[16A; 6 /k 2 C + 9kA 2 } 1/4 ■ exp[A v / ^C/(4/c 3 )] 



.(166) 



We can also consider relation f)l 30|) to the second order 
of approximation, which reads 



a(k) 



Ak 



IkC 2k 2 



(167) 



Constants c±(k) should, and can, still be chosen as in 
(fl65)l . in order to have \<j>\ ~ fc-3/2+(n 3 -i)/2 for large 

k. Once this new expression of a(k) is substituted in 
(I164p . it is instructive to plot the quantities \(f>(a(k) , k)\ 
and P = k 3 \(j>(k)\ 2 . In Figures ([S]), (TTUl) we can see these 
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FIG. 9: Field |0(fc)| versus k, in pre-inflation matter era (ar- 
bitrary units). 
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a ~ fc/V kC. Also here we may consider relation (| 129[) 
to the second (or further) order in k, 



i(k) 



2k 3 



(172) 



to get a better approximation. In Figures (fTT|) . (fT2|) we 
see plots of the field \<f>\ and of the primordial power spec- 
trum P(k). 
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FIG. 11: Field |</>(fe)| versus k, in pre-inflation radiation era 
(arbitrary units). 
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FIG. 10: Primordial power spectrum P(k) versus k, in pre- 
inflation matter era (arbitrary units). 



qualitative plots, where, for sake of simplicity, we arbi- 
trarily set the parameters A = B = C= k=1. 

For the pre-inflation radiation era, we start from 
Eq. (|154[) . With analogous steps we find 



1 



B 1 



and 



with 



K«) = ^ex P (---j 



v"{a) + F(a) 2 v(a) = 0, 



(168) 



F(a) 



9B 



B 



2 nl/2 



2 KCa 4 C*a 6 C 2 a 10 _ 
The WKB solution of Eq. (|154|) explicitly reads 



(170) 
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(169) FIG. 12: Primordial power spectrum P(k) versus k, in pre- 
inflation radiation era (arbitrary units). 



<f>k(a) = v k (a)p(a) = (171) 

e-"^ 4 v^C 2 i/a [ c+ (k)e lG ^ +c^{k)e- iG ^] 
[2nC 2 a 8 - k 2 Ca 6 + 9nBCa 4 + B 2 k} 1/a ■ exp[B/(4Ca 4 )] ' 

Again the functions c±(k) should be chosen such that 
the square bracket's content in the numerator goes as 
- fcC™*" 1 )/ 2 , so that |0| ~ fc-3/2+K-i)/2 for large 



Finally, in Figures ([13| , (|T4|) , we compare the plots for 
fields and power spectra, matter and radiation cases, in 
the same diagrams. 

We see that for both matter and radiation eras there 
is an exponential suppression of the low A:-modes. The 
matter diagram presents an interesting cusp, just be- 
fore dropping down, that is absent in the radiation di- 
agram. In the next section the effects of these features 
on the CMB spectrum will be further investigated, via 
deeper numerical analysis with the help of CMBFAST, a 
code specialized for cosmological simulations of the CMB 
power spectrum. 
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FIG. 13: Field |<^(fc)| versus k, for pre-inflation matter era 
(red/upper line), and radiation era (green/lower line). 
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FIG. 14: Primordial power spectrum P(k) versus k, for pre- 
inflation matter (red/upper line) and radiation (green/lower 
line) eras. 



E. Numerical computation of the primordial power 
spectrum 

1. Matter era 

For the numerical solution it is most convenient to use 
the equation of motion for the scalar field perturbation 
<f> written as a function of the variable a. So we can use 
equation (|131|) cast as 



[Aa + Ca 4 ]^ + 



r5 



A + 4Ca 3 



■ 4>k 



0. (173) 



The field perturbation 4>k{<i) is a function of k and a. 
In one of the previous sections we determined the param- 
eter of the matter contribution in the Friedmann equa- 
tion, A, to be 



j4 = 10- 1 /2^. 

v„ 



(174) 



After a radiation-dominance era and a subsequent pe- 
riod of matter-dominance before inflation, we can, with 
the help of Eq. (|103p and the assumption of a(t p ) = 1, 
compute the scale factor at the onset of inflation 



10 7 / 2 



When inflation starts, the two competing factors of mat- 
ter and inflation in the Friedmann equation (|75p must be 
of the same order of magnitude. This condition, or more 
rigorously, condition (|53")) , allows us to fix the coefficient 
C as 



C 



1 



2«L/ 



10 



'11 

v„ 



(175) 



For the numerical solution of the equation, we have to 
use the horizon crossing condition (|130l) to write a as a 
function of k. The equation for the scalar field pertur- 
bation, Eq. (|173j) . is evaluated numerically for a fixed ki, 
repeatedly for many different choices of ki, and the val- 
ues of <p(a(ki),ki) are assembled to form an evolution of 
the field perturbations in dependence of k. 
The boundary conditions for the numerical solution have 
been obtained in section IV Bl Using them to solve 
the differential equation (|173[) leads to the result shown 
in Fig. 1151 However, the power spectrum as a function 
of k is only a collection of data points. In order to be of 
any use for the CMBFAST code, it has to be given as an 
analytical function, which is obtained by fitting the data 
points with an opportune mathematical function. The 
function used for the fitting is 



P(k) = a 



f 



. Z..6 ' 

1 + — 



(176) 



where the parameters a, ...,g can be found as: 



parameter 


value 


a 


2.205 ■ 10~ 12 


b 


3.233 ■ 10~ 12 


c 


0.03 


d 


2.578 ■ 10~ 12 


e 


1.680 ■ 10" 9 


f 


1.593 ■ 10" 12 


g 


6.584 ■ 10" 14 



Fig. [15] shows the fitting function (full line) together 
with the numerical curve of the power spectrum (dashed 
line). 

With these results, it is then possible to continue with 
the evaluation of the CMB power spectrum by putting 
the fitting function into CMBFAST. 



2. Radiation era with totally evaporating black holes 

The equation of motion for the scalar field perturbation 
is Eq. (|132p . which can be cast in the form 



[B + Ca 4 ]^ 



2B 



- AC a' 



b' k + —fa = 0. (177) 
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FIG. 15: The numerical solution for the primordial power 
spectrum (dashed line) and its fitting function (full line) in 
the pre-inflation matter era scenario. 



The fitting in this case has been done using the func- 
tion 

P(k) = a + b- Arctan(c -k)+d- Arctan 2 (c • k), (182) 

with the parameters a,...,d to be adjusted by Mathe- 
matica. The resulting fitting parameters are given in the 
following table, and the curve (full line) together with 
the numerical result can be seen in Fig. 1161 



parameter 


value 


a 


-3.701 ■ 10 -18 


b 


1.261 ■ 1CT 1B 


c 


3.116 


d 


-5.935 ■ 10" 17 



The fc(a)-relation for the radiation case is given by 
Eq. (ll29[) . We shall compare the situation of the previous 
section, when massive remnants are left by the black hole 
evaporation, with the present situation, in which black 
holes are nucleated but then disappear again completely 
into radiation. 

Following the scenario of radiation in pre-inflation era 
developed in subsection IIV F[ we can determine the pa- 
rameters B and C in this case. When r spans the interval 
between the end of black hole evaporation, r r ~ 10 4 , and 
the onset of inflation, Ti n t ~ 10 8 , then the scale factor 
evolves in full radiation era as 



Pi lio-"] 



a(r) = a(t p ) ^— 

Reminding a(t p ) — 1 and r c ~ 
at the onset of inflation of 



1/6 



,1/2 



(178) 



10 3 we get a scale factor 



a(r m/ ) = 10 25 / 6 ~ 10 



4.16 



(179) 



The radiation content of the Universe is fixed at r r , the 
end of the black hole evaporation era. Therefore, since 
by then R H = 10%, 



B = p r {T r ) ■ a A (r r ) 
10 4 VCe P 



(180) 



(10 4 ) 3 U P 



[(10 3 )- 1 /6 (10 4 ) 2/3 ] 4 = lo2 /3^^ 



Here, as before, we choose again ( ~ 1. Thus, the infla- 
tionary parameter C can be fixed by requiring that it is 
of the same order as the radiation parameter B at the 
onset of inflation 



C 



B 



"inf 



10 



-16 



(181) 



Thus we now solve the equation for a scenario without 
GUP in the same way as before, and investigate the dif- 
ferences in the primordial power spectrum. Fig. ll6l shows 
the result (dashed line). 




FIG. 16: The numerical solution for the primordial power 
spectrum (dashed line) and its fitting function (full line) in 
the pre-inflation radiation era scenario, with completely evap- 
orating black holes (no GUP active). 



3. Radiation era without black holes 

The equation of motion for the scalar field perturba- 
tions can be directly taken over from the previous section 
- it describes the universe evolving from a state of radia- 
tion dominance into the inflationary period. Simply the 
parameters B and C in the equation have to be evaluated 
anew. 

The scale factor at the onset of inflation is, as previously 
mentioned, 



a(T inf ) = 10 7 / 2 



(183) 



At the Planck time, the radiation had Planck density, 
and so 



B = Pp( t p) ■ o- 4 ( t p) = 1 7T- 

Vp 



(184) 



22 



Thus, the inflationary parameter can be fixed as 

C=-^-~l(T 14 . (185) 

a inf 

Again the equation for a scenario without GUP is solved 
in the same way as before, and the resulting primordial 
power spectrum is shown in Fig.[T7] (dashed line). 

The fitting in this case is done using the same function 
as in the matter case, 

P(k) = a-- w + - ' (186) 

with the parameters a, g as in the following table: 



parameter 


value 


a 


-4.566 • 10" 15 


b 


3.383 ■ lCT 15 


c 


0.0089 


d 


-6.561 • 10" 15 


e 


152.074 


f 


3.783 ■ 10~ 16 


g 


3.564 • 10" 7 



Fit for the power spectrum 
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FIG. 17: The numerical solution for the primordial power 
spectrum (dashed line) and its fitting function (full line) in 
the pre-inflation radiation era scenario without black holes. 



VI. THE CMB POWER SPECTRUM 

In this section the previously calculated primordial 
power spectra, in the cases with GUP, without GUP, 
and without black holes, as well as the analytical re- 
sults for the approximated equations, will be fed into 
the CMBFAST code to obtain the CMB tempera- 
ture anisotropy spectrum that can be measured today. 
We will compare our results to the WMAP seven year 
data, and to the result for the CMB spectrum obtained 



by following the standard inflationary scenario without 
any pre-inflation era. 

The standard inflation (SI) model is in principle a very 
good fit to the data, considering all the different features 
it has to explain. Only the mode with I = 2 in the 
CMB spectrum is very low in comparison to the follow- 
ing data points. If this is not simply a statistical feature 
but the indicator of a new physical phenomenon, then 
the SI model has to be modified in order to satisfy the 
drop at the low I modes. For now, only the I — 2 mode 
can be used to construct models with a suppression at 
low I modes, but in the far future it will be possible to 
tell whether the tendency of the power spectrum to drop 
will continue further, whether it will stay at a lower but 
constant level or whether it might even rise up again. 

6000 
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FIG. 18: The binned CMB temperature anisotropy spectrum 
as measured by the WMAP satellite, in comparison with the 
standard inflationary scenario. 

To produce the CMB power spectrum corresponding to 
the scenario of standard inflation, we assumed a scalar 
spectral index of n s = 0.963 and a running of the index 
a n ^ = —0.022, i.e. the most recent result of the WMAP 
observations |45| . 

In the calculation of the CMB power spectrum, there is 
one parameter that can be varied, the number of e-folds 
of inflation. The total number of e-folds of inflation (from 
the start, when the mode k% left the horizon, to the end 
of inflation) can be given (see (46|) as 

Ntot = N(k )+ln^y (187) 

fco is the currently largest mode within the horizon, 

fc = 0.002 hMpc- 1 . (188) 

This is the pivot scale for the wavenumber that the 
WMAP team has been using in constraining inflation 
models from their data [i^, [43, EH . N(ko) is the number 
of e-folds from the time during inflation when this mode 
fco crossed outside the horizon, 

N(k Q ) = ln(A) • (189) 
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Observations can only constrain the number N(ko), as 
currently the mode k reenters the horizon, but no in- 
formation can be given about whether there were more 
c-folds AN of inflation before fco exited the horizon dur- 
ing inflation. The constraint on N(ko) stated in [46[ is 

7V(fc )=54±7. (190) 

Usually in the power spectrum fc is normalized to fco, 
and the power spectrum is taken at values P(^-), where 
fco = 0.002 hMpc . By dividing fc in the expression 
of the power spectrum by numbers smaller than fco, we 
can add more e-folds to inflation accounting for the time 
before fc exited the horizon: 

Ar -=Ki) + K|)=K|)- < 191 ) 

Instead of taking the power spectrum as before normal- 
ized over fco, we take it as Pf^-Y where any number 
with ki < fco is possible. 

So, varying fcj is equivalent to adding e-folds AN to the 
experimentally constrained number iV(fco) = 54 ± 7. 



A. Results from the numerically computed 
primordial power spectrum 

There are three cases to present from the numerical 
calculations. 

For a matter-dominated era before inflation the CMB 
power spectrum as obtained by CMBFAST can be seen 
in Fig.IHl 

For the case when the GUP is turned off, the CMB 
power spectrum can be seen in Fig. 1201 

For the case when there are no black holes at all, only 
pure radiation, the CMB power spectrum can be seen 
in Fig.HTJ 



B. Results from the analytical solutions for P(k) 

As regard the first and the second case of the list of 
the previous paragraph, the power spectrum obtained by 
WKB solving the approximated differential equation for 
a pre-inflationary matter and radiation era was fed into 
CMBFAST as well. The left picture in Fig.[22] shows 
the outcome for the CMB temperature anisotropy spec- 
trum in the case when the GUP is valid and thus the 
pre-inflation era is matter-dominated (we use the coeffi- 
cients A,C from Eqs. (|1741 I175p h whereas the right pic- 
ture in Fig. [22] shows the results in the case where GUP 
does not hold, which corresponds to a pre-inflation ra- 
diation era with completely evaporating black holes (the 
coefficients B, C from Eqs. (|180l I181[) have been used). 



C. Comparison 

To have a better impression of how the three scenarios 
compare to each other, in Fig.[23]there are several graphs 
showing three curves for different values of AN. Full 
lines represent the cases with GUP, dashed lines the cases 
without GUP, and the dashed-dotted lines the case of 
pure radiation, where no black holes existed at all. 

D. x 2 Calculations 

The goodness of the fit of each of the different curves 
can be quantitatively expressed by a x 2_v alue, calculated 
from the formula 

* = ATE C 2 1 . (192) 

i=i 1 

where Di is the i — th data point, Tj is the correspond- 
ing value calculated by the model, and Ci is the error 
bar of measurement for the i — th data point. The theo- 
retical model, calculated for different numbers of e-folds 
are compared with the WMAP 7- year-results for the un- 
binned CMB temperature spectrum. The error bar of the 
measurement is provided along with the WMAP spec- 
trum data. 

The % 2 -value for the SI model calculated according 
to Eg. (1192[) is Xsi = 1-154. For our model of matter- 
dominance in the pre-inflation era, the x 2_v alues for dif- 
ferent number of e-foldings are given in the following ta- 
ble. 



AiV 





0.693 


1.386 


1.792 


2.079 


2.996 


3.584 


4.094 


x 2 


4.801 


1.69 


1.206 


1.353 


1.463 


1.378 


1.171 


1.172 



For the scenario of radiation-dominance in the pre- 
inflation era, with totally evaporating black holes, the 
X 2 -values for different number of e-foldings are given as 
follows. 



AN 





0.693 


1.386 


1.792 


2.079 


2.996 


3.584 


4.094 


x 2 


1.121 


1.132 


1.141 


1.144 


1.147 


1.151 


1.152 


1.153 



For the scenario of radiation-dominance without the 
presence of black holes, the x 2 -values for different number 
of e-foldings are given as 



AN 





0.693 


1.386 


1.792 


2.079 


2.996 


3.584 


4.094 


x 2 


1.562 


1.33 


1.17 


1.157 


1.156 


1.155 


1.154 


1.154 



Fig. [M] should help to make the numbers more under- 
standable. It shows the x 2_va lues for the three models 
over the number of e-foldings AN added to the standard 
number of 54 e-folds. The horizontal line represents the 
value for the SI model. 
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FIG. 19: 

The CMB power spectrum for a pre-inflation matter era, for various cases of AN. Overall view and zoom into the lower 

multipole region. 




FIG. 20: 

The CMB power spectrum for a pre-inflation radiation era, for various cases of AN. Overall view and zoom into the lower 

multipole region. 
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FIG. 21: 

The CMB power spectrum for a pre-inflation radiation era without any black holes, for various cases of AN. Overall view and 

zoom into the lower multipole region. 
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FIG. 22: 

The CMB power spectrum from the analytic solution of the approximated differential equation with matter dominance and 

radiation dominance, for various cases of AiV. 






FIG. 23: 

Comparing the CMB power spectrum for pre-inflation matter (full lines), radiation with black holes (dashed lines) and 
radiation without black holes (dashed-dotted lines) eras for AiV = 0.693147, AiV = 2.07944, AiV = 2.99573, AiV = 3.58352 

and AiV = 4.09434. 
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FIG. 24: The distribution of x 2 -values for a pre-infiation mat- 
ter and radiation era, for various cases of AiV. Circles repre- 
sent the matter model, whereas squares stand for the radia- 
tion model and rotated squares denote the case with radiation 
only. The horizontal line denotes the value of the SI scenario. 



VII. CONCLUSIONS AND OUTLOOK 

In this work, we investigated the effects of an era be- 
fore inflation on the CMB power spectrum measured to- 
day. We utilized the phenomenon of black hole nucle- 
ation from quantum fluctuations of the metric in very 
early times to argue for the existence of several thousand 
micro black holes in the pre-inflation era, which cause the 
universe to be matter-dominated from about r ~ 10 3 t p , 
until the onset of inflation (t ~ 10 6 £ p ). By setting up 
the Friedmann equation of the universe evolving from 
matter dominance to an inflationary phase it is possible 
to calculate the power spectrum of primordial fluctua- 
tions of a scalar field $ living in this scenario, and then 
process this primordial power spectrum by CMBFAST 
code to yield the CMB temperature anisotropy spectrum 
measured today. As an alternative to the matter domi- 
nance scenario, we also investigated the implications of a 
radiation-dominated era before inflation by relinquishing 
the claim of the Generalized Uncertainty Principle, stat- 
ing that black holes can only evaporate down to Planck 
size. These two cases have been calculated both numer- 
ically and in analytical approximations. A third case 
is presented, in which black holes never existed and the 
pre-inflation era is purely radiation-dominated. 

From the overall analysis that has been done on the 
three different scenarios, only the model with matter- 
dominance in pre-inflation era is really successful in the 
suppression of the I = 2 mode, and incorporates the de- 
sired effect on the power spectrum very well. It asymp- 
totes to the standard inflationary model for higher AN, 
which is expected as with increasing duration of inflation 
the effects of a pre-inflationary era are shifted to larger 
scales and thus would be expected to influence the scales 
that are still to enter the horizon. 

For the case of radiation in pre-inflation era, with black 



holes evaporating until zero mass (no GUP), the result is 
rather poor - only one of the curves shows a suppression 
in the lower modes, all of them actually lie above the 
result given by a standard inflationary scenario without 
a pre-inflation era. The model shows a very good accor- 
dance for the high I region. 

In contrast to that, the model with only radiation in 
the pre-inflationary era turns out to be a little better; 
two of the curves show a suppression, while the others 
lie higher than the curve obtained by the standard infla- 
tionary scenario. This can be traced back to the shape of 
the primordial power spectrum, which in this case rather 
resembles the matter case than the radiation case with 
black holes. 

The CMB power spectrum, produced with the ana- 
lytical solutions obtained by approximating the differen- 
tial equation, looks very similar to the standard inflation 
picture with only slight deviations. But on the other 
hand, we know that the WKB analytical solution has 
been mainly used as a qualitative guide to choose the cor- 
rect numerical solution among several possibilities, and 
so we cannot expect from it a perfect fitting of the data. 
Especially because we pushed the analytical approxima- 
tion to the second order only. Of course, it is important 
to investigate the differential equation analytically as to 
obtain the correct boundary conditions for the numerical 
simulation and to firmly support the correctness of the 
numerical result. 

Although the suppression of the lower modes in the 
case of pre-inflation matter dominance is there, the over- 
all fit to the current data turns out to be not so successful. 
A drawback of the matter model definitely is the bad fit 
for large I when only a few e-folds are added; the behav- 
ior for large I becomes better with increasing AN, but 
it is not very good for small AA^. The reason why the 
model without GUP is so much better in the overall fit 
than the other two cases can be found in the fitting pro- 
cess for the primordial power spectrum. The shape of Pk 
in the matter model and in the pure radiation is more 
eccentric and harder to be fitted than in the case of radi- 
ation without GUP, where the power spectrum is rather 
smooth and the arctan-fit matches the curve quite well. 
In the matter case, the fitting function is more compli- 
cated and the quality of the fit is definitely worse. This 
of course influences the x 2 " va hie of the model, which is 
quite large for the smaller AN. A stronger suppression 
on several of the low I modes also leads to a larger devia- 
tion from the modes with 1^3. The SI model is without 
doubt the best fit in general; however, it doesn't capture 
the drop at the I — 2 mode. The model with matter in 
pre-inflation might not be the best fit on all scales, but 
in the future, with further modes on larger scales than 
I = 2 being suppressed, the model might become more 
successful than the standard inflation scenario. For sure, 
if the suppression of the I — 2 mode is to be continued 
with a suppression on even larger scales, the extension 
of standard inflation to having a pre-inflationary era is 
required, otherwise the drop for the lower modes will re- 
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main unexplained. 

The quality of the results has been put into numbers 
by a \ 2 analysis, which is the best way in the current 
situation to give a solid statement about the success of 
the three models today. 

In terms of this analysis, the pre-infiationary matter 
era is disfavored compared to the SI model. We could 
think that the roles might be really different in the far 
future, and the modes that are still to enter the hori- 
zon might be more successfully described by the mat- 
ter era scenario. However, even without looking too far 
into the future, but keeping our mind on the present, 
the pre-inflation matter era model seems to be the only 
one, among those here studied, able to capture and de- 
scribe the low I modes suppression. Further refinement 
of the model are obviously in order. And in principle, if 
a better matching with the observations will be achieved, 
the model can serve to directly check the validity of the 
Generalized Uncertainty Principle, as relation (|107[) ex- 
plicitly suggests, and more indirectly, of the Holographic 
Principle. These are surely intriguing avenues for future 
research. Only with time it will become clear whether 
the model with matter-dominance in pre-inflation era is 
superior to standard inflation in its success to explain the 
CMB power spectrum. 
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Appendix 1 

Here we compute the nucleation rate n* for micro black 
holes, i.e. the number of micro black holes of critical mass 
M created in a thermal bath of gravitons, via gravita- 
tional instabilities of hot flat space, per unit volume per 
unit time. Essentially, we follow the procedure detailed 
in Ref. [5j. Our discussion, as in that reference, will 
be based on the standard Heisenberg principle. As said, 
the GUP will be implemented in our argument only by 
considering the cutoff imposed on minimum masses and 
maximum temperatures. Critical mass and temperature 
are linked by the relation (fl4)) . which in standard units 
reads 



M 



he 3 



8TrGk B T 



(193) 



The probability for a quantum vacuum fluctuation 
to produce a black hole of critical mass M is 
exp(— AF/(k B T)), where AF is the change in the free 
energy of the system with T and V held fixed. Now 
AF = F — F g , where F is the free energy of the black 
hole and F g is the free energy of the thermal gravitons 
displaced by the black hole. F is related to the rest en- 
ergy of the black hole E = Mc 2 by 



E = F-T d -l 
dT 

This gives us a differential equation for F 



dF 
dT 



F 
T 



where, from (| 193[) 



E(T) = Mc 2 
Integrating F'(T) we find 



F{T) 



E(T) 
T 



he 5 
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(194) 



(195) 



(196) 



(197) 



To compute the thermal free energy F g of displaced gravi- 
tons we need, as it is clear from Eq. (|195l) . an expression 
for the total energy E of such thermal gravitons. This 
can be obtained from Eq. (f24|) . explicitly rewritten as 



E 9 TOT {V) 



15 c 3 h 3 



VT 4 



(198) 



where we chose the greybody factor for gravitons T g = 1, 
and we dropped the correction function A(/3,T), as the 
effects of GUP are considered only through the cutoff 
on masses and temperatures. Since the volume of the 
displaced gravitons coincides with that of the black hole 
(47i\Rf/3, with R s = 2GM/c 2 ), we have 



^TOT (T) 



k B T 
720 



Therefore Ea.([l95"| yields 



^(T) = -^logJ 



720 



(199) 



(200) 



So, finally, expressing things in Planck units we have 



AF 



1 



1 



k R T IGttO 2 720 



logO 



(201) 



Note that, in the range of interest, namely for < O < 1, 
we have | logO/720| <C l/(167r© 2 ), so the second term in 
(|201[) can be neglected. 

Knowing the probability of one statistical fluctuation, 
cxp(— AF/(k B T)), we have to estimate the density for 
such fluctuations. Consider the fluctuations on the small- 
est scale possible, namely with a wavelength A m ; n = a£ p , 
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a of order 1. Then, imagining a cubic lattice, the number 
of statistical fluctuations in the unit volume is 



n 



Amir 



(202) 



and the number of fluctuations per unit volume able to 
produce a black hole of critical mass M will be 



rig exp 



AF \ 



An 



exp 



1 



167t9 2 



(203) 



The number of micro holes with critical mass M, created 
per unit time per unit volume, can be therefore computed 
as 



dn* 
~~dT 



/ 1 



\8ir0 3 / dr 



(204) 



The value of d<d/dr can be obtained using Eqs.(fT4"l) and 
(|2"5)) . with the correction function A(/3, T) = 1, and T 7 = 
1. We have 



dQ 2tt 2 
rf7 _ 15 



e 4 



So finally 



dn* 7r 
^ = "*60 9 



(205) 



(206) 



In the statistical probability of one fluctuation, we should 
also include a term for the quantum correction of the free 
energy of the black hole. It can be shown [4{| that for a 
Schwarzschild metric we have 



j^quantuvn 



212 
'~45~ 



log 



[AC 



(207) 



where /i is a regulator mass of the order of the Planck 
mass. The final formula for the number of micro holes 
with critical mass M, created per unit time per unit vol- 
ume, reads, in Planck units, 



dn* 
~dT 



-eU- 

60 V28 



(208) 



A, 



7T_e /_a v 
60 \2e/ 



exp 



167t6' 



Since in Planck units i p = p, = 1 the previous formula 
can be usefully rewritten as 



dn* 
~dT 



87T 



15 • 64 7r 



<d 45 exp — 



1 



167r9 2 



(209) 



which agrees with the numerical pre- factor of Ref. [B[ since 
we chose a = it ■ 2~ 32 / 135 ~ 2.66. 



Appendix 2 

In this Appendix we construct the cosmic scale factor 
a as a function of the conformal time rj. Once we have 
a (n), we shall be able to write equation (|120p for Vk{rj). 

In order to arrive to the function a(rj), two different, 
but equivalent, procedures can be specified. 

I) When the scale factor a(t) is a known function of 
the cosmic time t, then we can compute (in principle, 
explicitly) the function n(t) 



n= I dn = 



dt 



v = v(t). 



Hence, inverting the last relation we get 

t = t(jf) a= a(t(n)) = a(n) . 



(210) 



(211) 



There exists however also another alternative procedure. 



II) Sometimes the integral in (|210l) is not easily doable, 
and moreover the object we are usually interested in is 
the function a(rj), and not rj(t). Such function can be 
directly computed by re-writing the equation of motion 
for a in terms of the conformal time rj, instead of the 
cosmic time t. The equation of motion (|72p reads, in 
conformal time, 



/ da 



k (Ca 4 + Aa + B) 



da 



^/Ca* + Aa + B 



i/Hdr) 



(212) 



(213) 



The function a(r]) can be obtained directly from the in- 
tegration of the previous equation. 

In the following, we shall integrate equation (1213[) in the 
two cases of our interest: pre-inflation radiation era, and 
pre-inflation matter era. 

Pre-inflation radiation era: In this case there is no 
matter, therefore A = 0. Equation (|2 13[) reads 



s/Kdr) = 
With the substitution 



da 



VCa* + B 



1/4 



the equation becomes 

(K 2 BC)^ 4 dn = 



dx 



We can now make use of the formula (see [50] ) 

dx 1 „ , 

= = - F(a,k) 



(214) 



(215) 



(216) 



(217) 
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where F(a, fc) is the elliptic integral of the first kind 

dfx 



F{a,k) , ^_ 

'o \/ 1 — k 2 sin 2 jj, 

In our specific case, Eg. (|217l) . we have 



(218) 



k = -r 



\-x 2 
TTx~ 2 



Then 



2(n 2 BC)^ V = r~= 

Jo J\ 



dfx 



(219) 



(220) 



yl — k 2 sin 2 \i 
Inverting the last integral, we get the Jacobi amplitude 



a = am[2 ( K 2 BC) 1/4 r) 
and, because of relation (|219j) . 
1-x 2 



(221) 



1 + x 2 



cosam[2 (k 2 BC) 1/4 t] 



2(k 2 BC) 1 / 4 t 1 , 



1 

71 



(222) 



where cn is the Jacobi cosine- amplitude (see again [501 ] 
for definitions and properties). Reminding the relation 
(|215p between a and x, finally we can write 



1 — cn 


2(k 2 BC) 1 ^7 1 , 


l 


1 + cn 


2(K 2 BC) 1 / i n 1 


1 

y/5. 



(223) 



It is interesting to check the small rj limit (or, equiva- 
lently, the small C limit) of the previous relation. Con- 
sidering the MacLaurin expansion for cn 



we have 



cn(u,fc) = 1- -u 2 + 0{u 4 ) 



B_Y 2(k 2 5C)V 



1 + 1 



(224) 



kBt] 2 (225) 



which is the well known expression for the scale factor 
a(rj) in pure radiation era. 

To build equation (|120[) for Vk we need to compute 
a" (rj) I a(rj) . It is a bit laborious, but however the result, 
for the parameter fc = 1/ \/2, is 



a[ri) 2 



1 — cn 



1 + cn 



0< 7l) ' 



(226) 



with P = 2{k 2 BC) 1 ' 4 . 

Therefore the equation for Vk reads 



4 



2 /3 2 l~cn(/fy) \ 

A — — • t-t-^- I Vt- = 



1 + cn(/3?7) / 



(227) 



where capital K is the cosmological perturbation wave 
number and has nothing to do with the elliptic functions 
parameter fc. This equation has only a resemblance with 
the Lame' equation, but unfortunately, differs from it in 
a fundamental way. 

Pre-inflation matter era: In this case we don't have 
radiation, i.e. B = 0, and the equation for the conformal 
scale factor a(rj) reads 



\/~Kdn ■ 
Using the substitution 



da 



VCa 4 + Aa 



a = x 



1/3 



we have 



(K 3 A 2 Cy/ 6 dr, 



dx 



^x{l+x3) 



(228) 



(229) 



(230) 



Again with the help of |5l| we can make use of the for- 
mula 



dx 



1 



^/x(l + x 3 ) </3 



(231) 



where F(a, fc) is the usual elliptic integral of the first kind 
(|2"Tg]> . but now 



y/2 + y/3 



a = arccos 



1 + (1 - x/3> 
1 + (1 + V3> 



Then 



^3 (k^C) 1 / 6 ^ 



dfi 



(232) 



(233) 



lo \/\ — k 2 sin 2 n 
and the Jacobi amplitude reads 

a = am[^3 (k 3 A 2 C) 1/6 rj\ (234) 
Reminding relations (|232l) and (|229[) . finally we have 

the solution a {rj) expressed in terms of Jacobi cosine- 
amplitude 



a{rj) 
where 



.4 



1 — cn[/3rj, fc] 



C) ( A /3-l) + (V3 + l)cn[/3r ? ,fc] 



(235) 



[3 = </3(K 3 A 2 Cy/ 6 r, ; fc = ^ 2 ± ^ . (236) 

The limit for small rj, or small C, can be computed with 
the help of Ea. (|224p . and reads 



airy) = ^ rj 2 



(237) 
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which is the known expression of the conformal scale fac- 
tor in pure matter era. 

As for the construction of the equation for Vk, we have, 
after a somehow long calculation, 



1 



V3/3 2 0(17) 



The equation for Vk can be henceforth explicitly writ- 
ten down, although it still results to be only "similar" 
to the Lame' equation, but not exactly of that known 
kind. Progress on the analytical exact solutions of such 
equations will be reported in future work. 



(238) 



(1 - V3)(V3cn(/37?) -!) + (! + V3)(V3 + cnQ^cu 2 ^) 
[(^-l) + (V3 + l)P(l-cn(^)) 
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